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CHAPTER  1 


INTRODUCTION 

During  the  past  decade,  substantial  investments  have  been  made  in  the  development 
of  compact  goo  discharge  pumped  CO2  lasers  for  military  and  civilian  applications.  The 
DOD  community  and  SDIO  in  particular  is  developing  compact,  lightweight  CO2  lasers 
for  airborne  and  space-based  radar  applications  directed  toward  target  ranging,  imaging, 
and  discrimination.  These  CO2  laser  systems  are  being  developed  in  the  LOWKATER 
and  LICD  programs. 

The  three  major  components  of  a  CO2  laser  radar  are  the  laser  head  and  flow  loop, 
the  pulsed  power  supply  and  the  optics  and  beam  control  system.  In  this  proposal  SRL 
will  present  a  novel  self-sustained  discharge  concept  that  should  result  in  the  stable  and 
efficient  extraction  of  large  specific  energies  (50  J /liter  atmosphere)  for  pulse  lengths  as 
long  as  100  ps.  This  factor  of  five  increase  in  the  energy  per  unit  volume  will  enable  laser 
heads  to  be  fabricated  which  are  factors  of  4  to  5  times  smaller  in  volume  and  weight  than 
with  conventional  discharge  approaches.  Further  the  smaller  cross  section  of  the  optical 
beam  will  result  in  a  significant  reduction  in  the  size  and  weight  of  the  optics  and  beam 
control  system. 

The  new  laser  discharge  concept  proposed  by  SRL  is  compatible  with  SDIO  applica¬ 
tions  which  require  efficient  operation  of  compact,  lightweight  CO2  lasers  at  kilowatt  power 
levels.  This  concept  is  based  on  a  new  discharge  method  which  will  insure  the  macroscopic 
and  microscopic  stability  of  the  discharge  at  high  pump  power  density  and  for  long  pulse 
lengths.  This  discharge  concept  promises  pulse  lengths  of  up  to  100  microseconds  du¬ 
ration.  scalability  to  multikilowatt  high  average  power,  high  volumetric  efficiency  (up  to 
50  joules/liter-atm),  and  high  electrical  efficiency.  This  new  discharge  method  utilizes 
spatially-uniform  x-ray  or  UV  preionization  of  the  laser  medium  followed  by  a  discharge 
pulse  which  initially  supplies  a  voltage  across  the  discharge  electrodes,  approximately  equal 
to  two  or  three  times  the  sustaining  voltage,  to  avalanche  the  electron  density  uniformly 

1-1 


SCIENCE  RESEARCH  LABORATORY 


to  1012  electrons/cm3.  The  energy  to  drive  the  discharge  is  supplied  by  an  all  solid  state 
rf-modulated  driver  with  high  incremental  source  impedance;  that  is,  the  driver  acts  as  a 
current  source  to  incremental  changes  in  impedance.  In  this  concept,  the  current  source 
drive  prevents  discharge  instabilities  which  increase  the  electron  density  and  therefore  the 
current  density  in  a  spatially-uniform  manner.  In  addition,  the  oscillating  discharge  voltage 
prevents  localized  electron  density  avalanche  by  controlling  the  ionization  and  recombina¬ 
tion  during  the  peaks  and  valleys  of  the  RF  waveform  and  by  limiting,  with  inductance, 
the  localization  of  discharge  current.  A  common  example  of  such  a  localized  instability  in 
a  flowing  CO2  laser  are  arcs  that  occur  in  the  downstream  side  of  the  laser  cavity  where 
the  laser  mixture  is  hotter  and  residual  ionization  is  generated  uniformly  at  the  electrode 
edges.  Present  methods  to  minimize  this  problem  rely  on  increasing  the  gas  flow  rate 
to  reduce  the  ionization  level  at  the  edges  of  the  discharge  electrodes  or  utilizing  electron 
beam  preionized  discharges.  Both  of  these  solutions  are  unacceptable  if  a  reliable  compact, 
lightweight  CO 2  laser  is  required.  Increasing  the  gas  flow  rate  will  increase  the  weight  and 
volume  of  the  laser  flow  loop.  Electron  beam  preionized  laser  discharges  require  a  thin  foil 
which  separates  the  vacuum  where  the  electron  beam  is  formed  from  the  high  pressure  laser 
mixture.  This  foil  can  create  reliability  problems  in  cases  where  long  duration  operation 
is  required  without  maintenance.  In  the  remainder  of  this  section,  existing  self  sustained 
techniques  for  discharge  pumping  CO 2  lasers  will  be  reviewed. 

The  x-ray  or  UV  preionized  discharge  method  utilizes  x-ray  or  UV  radiation  to  provide 
a  low  level  of  spatially  uniform  electron  density  which  can  subsequently  be  avalanched  uni¬ 
formly  by  applying  sufficiently  high  electric  fields.  This  avalanching  discharge  method  can 
typically  be  made  to  be  stable  for  long  single  pulse  operation  since  the  spatial  preionization 
contours  can  be  carefully  controlled  so  that  discharge  arcing  does  not  occur  at  the  edges 
of  the  electrodes.  However  repetitive  pulse  operation  requires  flowing  the  laser  medium 
which  carries  ionization  to  the  downstream  edges  of  the  discharge  electrodes.  The  electric 
field  enhancement  at  these  edges  can  then  trigger  discharge  instabilities.  Therefore,  this 
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discharge  concept  may  not  be  compatible  with  development  of  compact,  efficient  CO2  laser 
devices  for  space  based  applications. 

CO2  waveguide  lasers  have  proven  to  be  useful  and  compact  devices  for  various  tactical 
military  applications.  A  major  advantage  of  waveguide  devices  is  that  they  do  not  rely  on 
flowing  the  gas  to  remove  the  waste  heat.  The  narrow  waveguide  channel  allows  for  heat 
removal  by  conduction  through  the  sidewalls  of  the  waveguide  channel.  Unfort unately, 
waveguide  lasers  have  been  operated  at  low  pressure  (<  100  Torr)  and  so  their  utility 
in  laser  radar  has  been  severely  limited.  This  low  pressure  limitation  arises  because  the 
discharge  has  to  be  wall  stabilized.  In  fact,  the  DOD  has  made  a  considerable  investment 
to  alleviate  this  constraint  on  waveguide  laser  operation.  In  addition,  pulsed  operation  of 
50-100  watt  CO2  waveguide  laser  is  typically  possible  only  with  Q-switching  which  reduces 
efficiency  somewhat  and  also  reduces  reliability. 

Science  Research  Laboratory  has  developed  a  new  CO2  laser  discharge  concept  which 
circumvents  the  limitations  on  C02  laser  performance  discussed  above.  This  new  discharge 
concept  can  provide  stable  discharge  operation  for  pulse  durations  up  to  100  ^sec  while 
providing  pump  power  density  consistent  with  specific  laser  energy  densities  of  50  J /liter 
atmosphere.  The  key  technical  issue  addressed  by  this  new  compact  CO2  laser  discharge 
concept  is  maintaining  discharge  stability  for  pulse  durations  up  to  100  fisec.  In  this  new 
discharge  concept,  stability  is  maintained  by  two  distinct  processes.  First,  a  current  source 
is  chosen  to  drive  the  discharge  to  prevent  spatially  uniform  ionization  instabilities  from 
terminating  the  discharge.  With  current  source  drive,  if  the  electron  density  increases 
uniformly,  discharge  conductivity  increases  and  the  voltage  across  the  discharge  decreases. 
Since  for  CO2  laser  discharges  the  ionization  rate  decreases  rapidly  with  decreasing  E/N, 
the  ionization  rate  falls  and  the  discharge  is  stable.  Of  course,  stabilizing  the  total  discharge 
current  is  a  necessary,  but  not  sufficient,  condition  for  discharge  stability. 

Even  with  constant  current  drive,  instabilities  at  regions  of  high  electric  stress  can  lead 
to  localized  discharge  arcing  and  impedance  collapse.  To  stabilize  the  discharge  against 
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instabilities  having  a  small  spatial  scale  size,  oscillations  are  imposed  on  the  current  and 
voltage  pulses  to  prevent  constriction  of  the  discharge  current.  This  is  accomplished  by 
choosing  the  oscillating  frequency  to  be  sufficiently  rapid  to  allow  the  inductance  of  the 
constriction  to  suppress  current  flow  until  the  discharge  voltage  has  dropped  below  the 
sustaining  voltage  where  recombination  dominates  ionization.  Therefore,  the  oscillating 
voltage  waveform  prevents  localized  electron  density  avalanche  and  subsequent  arcing  while 
the  current  source  drive  prevents  instabilities  which  are  of  large  spatial  scale.  By  enhancing 
discharge  stability,  the  pulse  duration  and  laser  energy  extracted  per  liter  can  be  increased 
and  the  weight  and  volume  of  the  laser  head  and  the  flow  loop  can  be  decreased  by  factors 
of  4  to  5. 
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CHAPTER  2 


TECHNICAL  OBJECTIVES  OF  PHASE  I 

Science  Research  Laboratory  has  successfully  met  all  the  objectives  of  Phase  I.  The 
primary  objective  of  this  phase  was  to  develop  a  two  dimensional  discharge  model.  This 
model  which  will  be  discussed  later  in  the  report  predicts  stable  operation  for  pulse  lengths 
as  long  as  100  /zs.  The  second  objective  of  the  Phase  I  effort  was  to  develop  a  conceptual 
design  for  a  repetitively  pulsed  CO2  laser  experiment.  SRL  has  completed  the  design  of 
a  laser  head  and  pulse  power  system  that  is  capable  of  operating  up  to  a  pulse  repetition 
rate  of  50  Hz.  The  details  of  this  design  is  also  presented  in  this  report.  The  SRL  two- 
dimensional  model  has,  of  course,  been  used  to  model  the  discharge  in  the  laser  head 
design. 

During  Phase  I  of  the  proposed  program,  a  comprehensive,  interactive  model  for  the 
new  CO2  laser  discharge  concept  was  developed  and  used  to  determine  regions  of  stable 
discharge  operation  and  efficient  C02  laser  performance.  In  the  model  it  was  assumed  that 
x-ray  or  UV  preionization  provides  a  spatiaily-uniform  electron  density  of  approximately 
106  -  108  electrons/cm3  which  is  then  increased  to  approximately  1012  electrons/cm3  in 
200  -  300  nsec  by  avalanching  during  the  front  end  of  the  main  discharge  pulse.  The 
main  discharge  pulse  maintains  the  electron  density  at  approximately  1012  electrons/cm3 
at  discharge  electric  fields  (and  therefore  electron  temperatures)  which  efficiently  pump  the 
CO2  laser  mixture.  The  discharge  voltage  is  assumed  to  be  sinusoidal  at  a  given  frequency. 

The  period  of  oscillation  is  chosen  to  be  short  compared  to  the  CO2  laser  transition  life¬ 
times  and  the  time  scale  on  which  ionization  instabilities  develop.  On  these  grounds,  it  is 
estimated  that  the  oscillation  frequency  should  be  between  10  and  102  megahertz.  The  up¬ 
per  limit  on  oscillation  frequency  is  constrained  by  the  requirement  that  the  discharge  skin 
depth  be  greater  than  the  transverse  dimension  of  the  discharge.  These  conditions  ensure 
a  continuous,  near  constant  laser  output  during  the  pulse  as  well  as  instability  suppression. 
An  analysis  which  treates  volumetric  (large  scale)  discharge  stability  will  be  presented  in 
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the  following  chapters.  The  remainder  of  this  chapter  discusses  the  consideration  that  goes 
into  choosing  the  RF  frequency. 

RF  modulation  of  the  discharge  pulse  will  have  an  impact  on  discharge  and  laser 
parameters  such  as  electron  density,  temperature  and  small  signal  gain.  These  effects  will 
be  discussed  briefly  below. 

The  electron,  ne,  production  and  loss  rate  is  given  by 

^  ~  S  +  V{Tit  -  cm2  -  0nt  (2.1) 

where  i /,•  is  the  electron  impact  ionization  rate,  a  is  the  electron  ion  recombination- rate  and 
0  is  the  attachment  rate.  For  a  self  sustained  discharge  the  source  function,  S=0.  Typically, 
CO2  laser  discharges  are  recombination  dominated.  Because  of  the  RF  modulation,  the 
discharge  electric  field  will  vary  in  time.  The  electron  temperature  and  recombination 
rate  will  also  vary.  As  the  electric  field  decreases,  the  average  electron  energy  decreases 
and  recombination  rate  increases  which  will  have  a  stabilizing  effect  and  allow  slightly 
constricted  regions  in  the  discharge  to  recover  to  the  ambient  electron  density. 

It  is  desirable  to  keep  the  electron  density  constant  in  time.  From  Eq.  (2.1)  it  is  clear 
that  for  0  <<  orne,  the  electron  density  will  be  constant  if  f  >>  ane.  Typically,  ne  ~  1012 
cm-3  and  ck  ~  10~6  cm3/sec.  So  if  f  >>  1  MHz,  ne  will  be  essentially  constant. 

The  RF  electric  field  can  also  lead  to  variations  in  the  average  electron  energy.  Reduc¬ 
tions  in  the  electron  temperatare  durng  minima  in  the  applied  RF  modulation  to  the  main 
discharge  pulse  allow  localized  regions  of  increased  electron  density  to  relax  during  the  RF 
pulse  to  insure  discharge  stability.  The  electron  temperature  will  vary  appropriately  if 
the  RF  modulation  frequency  is  /  <<  ktNa  where  ke  is  the  electron  impact  excitation 
rate  of  ground  state  of  N2  and  Na  is  the  N2  density.  Typically  10_9cm3/sec  and 

Na  ~  5  x  1018cm-3  and  therefore  f  <  2  x  10 10  Hz.  Since  the  RF  modulation  frequency  will 
certainly  be  lower  than  1010  Hz,  the  average  electron  energy  will  follow  the  modulation 
frequency  in  time  and  therefore  discharge  stability  will  be  enhanced. 
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For  the  small  signal  gain  to  be  constant,  f  >>  ke  Nc,  where  ke  (~  10~n  cm* /sec)1,1) 
is  the  vibrational  transfer  rate  from  N2  to  CO2  and  Nc  1018  cm_3is  the  CO2  number 
density.  So  for  a  constant  small  signal  gain,  f  must  be  greater  than  10  MHz.  The  streamer 
stability  criterion  is  consistent  with  this  modulation  frequency  so  the  small  signal  gain  will 
be  constant  which  is  important  for  maintaining  high  laser  power  extraction  efficiency. 

For  efficient  coupling  of  power  into  the  laser  discharge  it  is  important  that  the  impedance 
of  the  discharge  be  roughly  constant  during  the  RF  modulation  period.  The  power  reflec¬ 
tion  R  from  a  load  having  an  impedance  is  given  by  R  =  (l  -  Z0/Zi)2/(1  +  Z0/Zi)2 
where  Z0  is  the  impedance  of  the  RF  power  supply  and  Zi  oc  E/enevo  where  Z  is  the 
applied  electric  field.  From  a  Boltzmann  code  analysis  of  typical  C02  laser  mixtures  the 
drift  velocity,  vq  varies  approximately  as  (E)1-/2.  From  Eq.  (2.1),  ne  ~  uja.  i/j  is  a  very 
strong  function  of  the  electric  field  while  a  decreases  with  increasing  electric  field.  So  ne 
could  vary  quite  strongly  with  the  electric  field,  i.e,  ne  ~  E7.  Such  a  strong  variation  could 
result  in  power  reflection  from  the  discharge  and  the  overall  efficiency  of  the  discharge 
could  be  adversely  effected.  It  is  in  part  for  this  reason  that  the  RF  frequency  be  chosen 
such  that  the  electron  density  is  constant,  i.e.  f  >>  ane  ~  106  Hz. 

Finally  the  cross  sectional  area  of  the  discharge  will  be  constrained  by  the  skin  depth 
of  the  RF  portion  of  the  discharge  current.  For  a  given  RF  frequency  the  cross  sectional 
dimensions,  h,  must  be  smaller  than  a  skin  depth  if  the  current  through  the  cross  section 
is  to  be  uniform.  This  unequality  is  given  by 


where  cr,  the  conductivity  of  CO2  laser  discharges,  is  typically  2-3  x  10"2  mhos/M.  This 
constraint  has  been  discussed  quantitatively  in  conjunction  with  the  development  of  the 
criterion  for  streamer  suppression. 

From  the  preceding  discussions,  a  scaling  map  for  C02  can  be  drawn  as  shown  in  Fig. 

2.1.  This  figure  is  a  plot  of  RF  frequency  vs  cross  sectional  area  of  the  discharge.  As 
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discussed  earlier  the  RF  frequency  should  be  larger  than  30  MHz  for  the  small  signal  gain 
and  electron  density  to  be  constant.  On  the  high  end.  the  frequency  is  constrained  by 
the  skin  effect.  Also  shown  in  Fig.  2.1  is  the  expected  10.6  /i  laser  energy  that  can  be 
extracted  per  pulse.  From  this  figure  it  is  clear  that  this  pumping  technique  can  provide 
enough  laser  energy  for  laser  radar  applications. 


I 
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Figure  2. 1 :  Clearly  operation  at  modulation  frequencies  as  high  as  possible 
consistent  with  the  discharge  skin  depth  constraint  is  optimum 
for  early  suppression. 
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CHAPTER  3 


REVIEW  OF  THE  ONE  DIMENSIONAL  DC  VOLUMETRIC 

DISCHARGE  STABILITY 

3.1  INTRODUCTION 

The  stability  of  high  pressure  discharges  has  been  a  subject  of  major  concerns  in  the 
high  power  laser  community  over  the  last  fifteen  years.  Most  of  the  discharge  stability 
research,  especially  for  CO2  laser  discharges  concentrated  on  thermal  instabilities*1)  and 
their  effect  on  the  constriction  and  arcing  of  the  discharge.*2)  Vibrational  stabilization 
methods  have  been  proposed  to  suppress  these  instabilities.*3) 

The  thermal  instabilities,  by  their  nature,  can  only  develop  on  the  acoustic  (millisec¬ 
ond)  time  scales.  However,  in  many  cases*4)  the  arcing  occurs  on  a  much  faster  (micro¬ 
second  and  submicrosecond)  scale.  The  phenomenon  is  characteristic  of  systems  where 
the  electron  loss  by  attachment  is  present  and  the  two-step  electron  impact  ionization  is 
dominant*4).  Such  a  model  successfully  explained  the  experimentally  observed  temporally 
rapid  instability  scales  in  the  KRF  laser  discharges.*4) 

In  addition  to  the  KRF  discharge  case,  subacoustic  arcing  time  scales  have  been  also 
observed  in  the  high  power  CO2  laser  experiments.  An  example  of  such  arcing  is  in  con¬ 
nection  with  the  gas  recirculation  in  the  laser,  resulting  in  the  buildup  of  the  attaching 
oxygen  component  in  the  laser  gas  mixture.*5)  These  observations  could  also  be  explained 
within  the  two-  step  ionization  model.*4)  Although  the  details  of  the  kinetics  in  a  high  pres¬ 
sure  molecular  discharge,  are  complex,  the  concept  of  two-step  ionization  as  the  dominant 
ionization  source  in  the  discharge,  is  plausible.  The  basic  reason  stems  from  the  fact  that 
the  electron  temperature  in  these  discharges  are  relatively  low  (typically  1-2  electron  volts) 
and  very  few  electrons  in  the  tail  of  the  distribution  function  can  directly  ionize  the  ground 
state,  while  the  ionization  of  a  variety  of  available  excited  and  metastable  states  requires 
much  lower  energies  and  is  accessible  by  the  bulk  of  the  electron  energy  distribution.  It  is 
precisely  this  reason,  though  speculative  in  its  details,  but  almost  certain  in  reality,  which 
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makes  the  further  development  of  the  multi-step  ionization  discharge  theory  necessary. 

3.2  THE  STEADY  STATE  SOLUTION 


Consider  a  system,  shown  schematically  in  Fig.  3.1  with  a  DC  voltage  source  V0.  The 
voltage  is  distributed  between  the  external  resistor  R  and  the  plasma,  so  that 


v;  —  RI  +  LE 


(3.1) 


where  I  is  the  discharge  current,  L  is  the  plasma  length  and  E  is  the  electric  field  in  the 
plasma,  which  is  assumed  to  be  uniform  at  all  times.  In  writing  Eq.(3.1)  it  is  assumed  that 
the  cathode  fall  voltage  is  negligible  compared  to  LE  (a  situation  typical  of  high  pressure 
discharges). 

Three  regimes  of  operation  will  be  investigated: 

A.  Stiff  Voltage  Case  -  Here  R  ~  0  and  thus  the  field  in  the  discharge 


E  =  —r  =  Const 
Li 


(3.2) 


is  not  affected  by  a  possible  time  variation  of  discharge  characteristics,  such  as  the 
electron  density,  temperature  etc. 

B.  Stiff  Current  Case  -  In  this  case  R  — ►  oo  and  V0  — *  oo,  so  that  the  current 


,  V'o 

i  ~  =  const 

R 


(3.3) 


is  time  independent  and  fully  described  by  the  external  circuit. 

C.  Finite  External  Impedance  Case  -  In  this  intermediate  case  (0  <  f?  <  oc) 
both  the  current 

V0~LE(i) 


/(*)  = 


R 


(3.4) 


and  the  electric  field  E(t)  depend  on  the  details  of  the  time  evolution  of  the  dis¬ 
charge. 

At  this  point  it  will  be  informative  to  discuss  the  kinetics  of  the  self-  sustained  dis¬ 
charge.  The  densities  of  various  species  in  the  discharge  will  be  assumed  to  be  spatially 

3-2 

SCIENCE  RESEARCH  LABORATORY 


uniform,  and  the  discharge  is  described  by  the  following  system  of  rate  equations 

dn  £ 

dt 
dxi  m 


=  vnmne  —  an 2  —  ftne 


nr 


—  7neno  ~ 

at  r 


(3.5) 

(3.6) 


Here  ne  is  the  electron  density,  governed  by  the  ionization  of  an  excited  or  metastable 
state  having  a  density  nm  and  by  recombination  and  attachment  losses  (u,  a  and  ft  are  the 
corresponding  rate  constants).  The  density  of  the  excited  or  metastable  state  is  defined 
by  the  balance  between  the  processes  of  the  excitation  by  electron  impact  of  atoms  from 
the  ground  state  (of  density  na)  and  the  collisional  quenching  by  the  neutral  gas  (with 
characteristic  time  constant  r).  The  time  scales  r,(ane )  and  ft~l  axe  assumed  to  be  much 
shorter  than  the  acoustic  times  and  thus  the  neutral  gas  temperature  and  density  are 
constants. 

Equations  (3.5)  and  (3.6)  yield  the  steady  state  in  which  all  the  discharge  characteristics 
are  constant  and 


^  m  —  nmo  —  7o^’^eo^a 

fto 


ne  —  nt0  — 


(3.7) 

{voloTna  -  a0) 

The  question  of  whether  this  steady  state  can  be  reached  for  a  given  set  of  parameters  of 
the  external  circuit  (V'0,  R  and  L )  will  be  considered  next.  The  answer  to  this  question 
depends  on  the  regime  of  operation: 

A.  Stiff  Voltage  Case  -  Here  E  =  E0  is  given,  which  defines  all  the  rates  in  Eq. 
(3.S)  so  that  ne0  is  known.  This,  in  turn,  defines  the  current  density  j  —  ent0vi0 
in  the  plasma  ( Vd0  being  the  drift  velocity  of  the  electrons).  Thus  the  steady  state 
exists,  provided  of  course  (see  Eq.  (3.8))  the  quantity 

Q{E0)  =  i/„7 0rna  -  a0  (3.9) 

is  positive.  This  condition  restricts  the  accessibility  of  the  steady  state  and  needs 
a  further  discussion.  Consider  the  dependence  of  Q  on  E0 •  Typically,  the  ex¬ 
citation  and  ionization  rates  (7 oi^o)  are  increasing  functions  of  the  electric  field, 
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while  the  attachment  rate  a0  normally  decreases  with  E0.  Furthermore,  ^’ZqQ  — 
—0c(Eo— o)  <  0  and  with  the  increase  of  E0,  Q  —*  u0j0Tna  >  0.  Therefore  there 
exists  a  single  value  of  the  electric  field  E0  =  E *  at  which  Q(E')  =  0.  The  steady 
state  operation  of  the  discharge  under  the  stiff  voltage  conditions  ( E0  is  given)  is 
thus  only  possible  when 

E0  -  y-  >  E'  (3.10) 

Li 

B.  Stiff  Current  Case  -  Here  the  current  I0  is  defined  by  the  external  circuit.  On 
the  other  hand  from  Eq.  (3.S) 

I0  =  evioneoS  =  ~^T  =  W)  (3.H) 

where  5  is  the  cross  section  area  of  the  discharge.  Equation  (3.11)  now  defines  the 
steady  state  electric  field  value  E0.  A  Graphical  illustration  of  the  solution  of  this 
equation  is  given  in  Fig.  3.2,  where  the  qualitative  dependence  of  function  F  on 
E  is  shown  based  on  our  earlier  discussion  of  function  Q(E).  A  possible  minimum 
of  F  at  E  =  Em  is  also  shown.  Such  a  minimum  may  exist  at  high  values  of  E0  at 
whicu  the  ionization  and  excitation  rates  grow  slower  than  the  drift  velocity  vd  as 
the  electric  field  increases.  It  can  be  seen  in  Fig.  3.2  that  the  steady  state  exists  for 
any  current  provided  it  is  larger  than  Im .  It  is  noteworthy  that  the  operation 
near  E*  requires  rather  high  currents. 

C.  Finite  Impedance  Case  -  In  this  intermediate  case  Eq.(3.11)  is  replaced  by  (see 
Eq.  (3.4)  for  the  current) 

—  ~£E°  =  F(E.)  (3.12) 

The  graphical  solution  of  this  equation  for  E0  is  shown  in  Fig.  3.3.  One  observes 
that  for  given  values  of  V0  there  exists  a  critical  impedance  Rcr  such  that  when 
R  <  Rcr  two  steady  state  solutions  (E0 1  and  fJo2)  exist  for  the  electric  field,  while 
for  R  >  Rcr  no  steady  state  exists  at  all.  For  R  —  Rcr  there  is  a  single  steady  state. 
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F(E) 


.V  /  R 
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R  <  R 


■—  Vo/L 


Figure  3.3:  Graphical  solution  for  the  electric  field  Eo  in  the  discharge  under  finite  impedance 
conditions.  The  straight  lines  describe  the  function  (Vo  -  LEo)/R  for  different 
values  of  the  external  impedance  R. 
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The  value  of  Rcr  can  be  found  from 


L  =  -RctF'{Ect) 

where  F'  —  dF/dE0  and  Ecr  is  given  implicitly  by 

F  K  _  F(Ecr) 

CT  L  F'(Ect) 


(3.13) 


(3.14) 


This  completes  the  analysis  of  the  possiblity  of  a  steady  state  solution  at  various 
external  circuit  conditions. 

3.3  VOLUMETRIC  STABILITY  ANALYSIS 

At  this  stage  one  can  analyse  the  problem  of  the  stability  of  the  steady  state  discharge 
with  respect  to  a  small  uniform  (volumetric)  perturbations  of  the  electron  and  excited 
state  densities.  Accordingly,  one  looks  for  solutions  of  Eqs.  (3.5)  and  (3.6)  in  the  form 

ne(t)  =  ne0  +  8ne(t) ;  nm(t)  =  nmo  +  6nm(t)  (3.1) 

where  neo  and  nmo  are  the  steady  state  solutions  described  in  the  previous  section  and 
8nt.8nm  are  small  perturbations.  Again  consider  the  stiff  voltage,  stiff  current  and  finite 
impedance  cases  will  be  considered  separately. 

A.  Stiff  Voltage  Case  -  Since  the  electric  field  in  the  discharge  in  this  case  is  a 
constant  of  time  ( E0 )  all  the  rates  in  Eqs.  (3.5)  and  (3.6)  are  also  constants  of  time 
and,  after  the  linearization,  one  can  write 
d(8nt) 


dt 

d(8n,n) 

dt 


=  vnt08nm  +  (i mmo  -  2 an, 

c  8nm 
—  7naonc 


f3)8nt 


(3.2) 

(3.3) 


Assuming  solutions  8nm,8ne 
equation 


exp(  —  icot),  one  then  obtains  the  characteristic 


^  +  iu:(oine0  H — )  H —  =  0 

T  T 


The  roots  of  Eq.  (3.4)  are 


^’1,2  =  —^ 


G^eo  +  ~  I  ± 
r 
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and,  since  one  of  these  roots  has  Imu  >  0,  the  solution  of  Eqs.  (3.2)  and  (3.3) 
is  generally  unstable  [except  when  (3  —  0,  in  which  case  one  of  the  roots,  w2  = 
0,  indicates  a  marginal  stability].  It  can  be  also  seen  that  if  r  is  the  shortest 
characteristic  time  in  the  problem,  Eq.  (3.3)  yields  uj2  —  i f3,  i.e.,  the  growth  rate 
is  approximately  the  rate  of  attachment. 

B.  Stiff  Current  Case  -  Here,  the  current  I  =  ev<intS  in  the  discharge  remains 
constant  in  time,  regardless  the  variation  in  the  electron  density.  Therefore, 

Svj,nt0  +  SritVdo  =  0 


and,  since 


c  dvdo  c 

=  1e.6E 


the  electron  density  perturbation  results  in  the  perturbation  of  the  electric  field, 
given  by 

SE  a=  (3.6) 


v'jon«> 

where  Vdo  =  dvdo/dE0.  Typically  Vdo  >  0  and  thus  the  variation  of  the  electron 
density  leads  to  the  opposite  change  in  the  electric  field  E0  in  the  discharge.  This,  of 
course,  plays  a  stabilizing  role  via  the  appropriate  change  in  the  electron  production 
rate.  The  phenomenon  can  be  demonstrated  analytically.  Indeed,  in  the  case  of 
interest,  perturbed  Eqs.  (3.2)  and  (3.3)  should  now  be  replaced  by 


d(6ne) 

dt 


“b  /?0)<^77e 


+  Knmoneo  -  ct'0n\g  -  ld'0nt0)6E 


d(6nm)  ,  i  r  r, 

- r —  =  7 onaonf - 4-  7 0nant0dE 

at  r 


(3.7) 


(3.S) 


where  <f'  =  d£/dEa  for  £  =  vo,ao,  80  and  7„.  The  substitution  of  expression  for  6E 
transforms  Eqs.  (3.7)  and  (3.8)  into  a  system  of  linear  first  order  differential  equations  for 
6nt  and  6nm.  The  characteristic  equation  for  this  system  is 


u?  -f-  iujA  +  - - 2?  =  0 

r 


(3.9) 
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where 


^  u  ^ 

A  =  OcTleo  -f*  I  J~(^mo^o  n<0<20  0o) 


T  V, 


(3.10) 


Vdo  L  _  /  .  7o\  «>e<,  + 


B  =  -P  a0nav0nt0  -£  +  _ 


v0  To 


(3.11) 


The  roots  of  Eq.  (3.9)  are 


aq,2  =  -i  A±Ja*+<i(^-b' 


(3.12) 


Finally,  by  inspecting  Eqs.  (3.5  and  (3.11))  one  observes  that 

0  B  Q(Eo)F'(E0 ) 

t  eSrv'j. 


(3.13) 


and  therefore  fl/r  —  B  <  0,  if  F\E0)  <  0.  In  this  case  Eq.  (3.12)  yields  characteristic 
roots  with  Imu  <  0  and  thus  the  discharge  is  stable  in  the  decreasing  region  of  function 
F(E0)  (E*  <  E0  <  Em  see  Fig.  3.2.). 

C.  Finite  Impedance  Case  -  In  this  intermediate  case  (see  Eq.  (3.12)) 


—  =  ev<ineS  +  -77-  =  const(t) 
fC  K 


(3.14) 


Therefore,  instead  of  (3.6),  one  has 


c  r  _  _  _ _ 

d°Wdo  +  L/eRSnto)nt 


(3.15) 


This  shows  that  Eqs.  (3.9)  through  (3.12)  also  holds  in  the  finite  impedance  case 
if,  in  these  equations,  one  replaces 


do  *  V deff  —  *-’(io  "b  1/ eRSrit 


(3.16) 


As  R  — *  00,  Vd'fj  — ►  v'do  and  the  stiff  current  discharge  is  recovered  (except  for  the 
fact  that  V0  is  finite)  (Case  B).  If,  in  contrast,  R  — ♦  0  then  v'de^j  — ♦  00  and  B  — ►  0, 
which  is  the  unstable  case,  corresponding  to  the  stiff  voltage  conditions. 


SCIENCE  RESEARCH  LABORATORY 


3.4  NUMERICAL  EXAMPLE 


In  this  section  the  discharge  model  developed  so  fax  will  be  applied  to  the  specific  case  of 
a  CO2  laser  discharge.  Before  this  can  be  done,  however,  one  has  to  know  the  relevant  rate 
constants  for  the  COo  mixture  as  a  function  of  the  electric  field.  A  one  atmosphere  toted 
pressure,  3/2/1,  He/N2/C02  mixture  will  be  chosen  as  an  example.  The  rate  constants 
for  the  ionization  and  excitation  of  the  electronic  states  of  iV2  in  this  mixture,  the  fraction 
of  discharge  energy  that  goes  into  vibrational  excitation  and  the  drift  velocity  as  predicted 
by  the  Boltzman  code^  are  shown  in  Figs.  3.4  and  3.5.  Form  Fig.  3.4,  which  is  a  plot 
of  Te,  vj  and  the  efficiency  of  exciting  the  vibrational  levels  of  JV2.  It  is  clear  that  efficient 
CO2  laser  operation  occurs  between  5-10  kV/cm. 

Figure  3.5  shows  the  variation  of  v  and  7  (see  Eqs.  (3.5  and  (3.6))  as  functions  of 
the  electric  field.  Also  shown  in  Fig.  3.5  is  the  curve  for  the  ionization  rate  from  the 
ground  state  (17,).  The  curve  for  7  shown  in  Fig.  3.5  is  the  total  excitation  rate  constant 
for  all  the  electronic  levels.  This  was  chosen  since  it  is  difficult  to  identify  the  relevant 
metastable  or  electronic  state  that  will  be  subsequently  ionized.  It  is  probable  that  this 
state  will  be  a  high  lying  state  such  as  the  C  state  of  N 2.  The  ionization  rate  constant  u 
of  the  excited  electronic  state  has  the  same  shape  as  the  electron  impact  ionization  rate  of 
N2  (a  ''  its  magnitude,  however,  was  increased  by  102,  since  this  state  is  about  9  eV 
below  the  ionization  level  of  JV2  as  opposed  to  the  C  level  which  is  less  than  4  eV  below 
the  ionization  level.  From  Fig.  3.5  it  is  clear  that  va  «  v  in  the  region  of  interest  and  so 
for  any  reasonable  metastable  density  ( nm  >  1012  cm'3)  the  metastable  ionization  will  be 
dominant.  Fig.  3.6  shows  the  recombination  coefficient  as  a  function  of  the  electric  field. 
Measurements  of  the  recombination  coefficient  for  a  3/2/1  laser  mixture  has  been  made 
for  electric  field  strengths  between  2.5  and  5  kV/cm/6*  The  curve  shown  in  Fig.  3.6  makes 
use  of  these  measurements  and  assumes  an  almost  constant  recombination  rate  for  10  < 

E  <  15  kV/cm. 

3-11 

SCIENCE  RESEARCH  LABORATORY 


SCIENCE  RESEARCH  LABORATORY 


10' 


a 

cm^/sec 


10' 


10' 


5 


Figure  3.6:  The  dependence  of  the  recombination  rate  on  the  electric  field  in  a  1  atm, 
3/2/1  CO2  laser  mixture  discharge. 
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At  this  stage  the  stability  of  the  discharge  at  various  operating  conditions  .will  be 
analysed.  The  analysis  will  be  started  by  presenting  the  dependence  of  function  F(E)/S 
on  the  electric  field  (see  Eq.  (3.11)).  This  dependence,  is  shown  in  Fig.  3.7  for  the  case  of 
interest.  One  observes  that  E *  ~  S.05  kV/cm,  so  that  the  steady  state  discharge  operation 
is  possible  for  8.05  kV/cm  <  E  <  15  kV/cm  (see  the  discussion  in  Section  2  and  Figs. 

3.2  and  3.3).  The  graphical  solution  for  the  electric  field  Ea  in  the  steady  state  discharge 
is  also  shown  under  the  following  conditions:  V0  =  9  kV,  RS  =  400  Qcm2,  L  =  1  cm, 
r  =  1  n sec  and  0O  =  105  sec-1.  The  steady  state  electric  field  is  S.45  kV/cm  which,  by 
using  Eq.  (3.8),  and  data  in  Figs.  3.4  through  3.6,  yields  the  steady  state  plasma  density 
neo  =  1.07  •  1012cm-3.  Prior  to  studying  the  stability  of  the  discharge  at  these  conditions, 
the  stiff  voltage  situation  is  considered.  The  analysis  in  Section  3  shows  that  the  discharge 
is  unstable  in  the  stiff  voltage  case  for  fi0  ^  0.  This  unstable  behavior  is  demonstrated  in 
Figs.  3.8  and  3.9,  where  the  numerical  solutions  of  Eqs.  (3.5)  and  (3.6)  are  shown  for  E0 
=  8.45  kV/cm  =  Const(t)  (stiff  voltage  condition)  and  two  different  initial  values  for  the 
electron  density.  The  electric  field  E0  defines  various  rates  in  Eqs.  (3.5)  and  (3.6)  which 
remain  constants  in  the  calculations,  consistent  with  the  stiff  voltage  conditions.  Figure 
3.8  shows  the  case  when  ne(0)  =  nm(0)  =  5.10ncm-3.  This  initial  density  is  lower  than  the 
above  mentioned  steady  state  density  (neo  =  1.07  •  1012cm-3).  Figure  3.8  shows  that  the 
electron  density  decays  in  time,  in  this  case,  and  the  discharge  is  quenched.  Figure  3.9,  in 
contrast,  shows  the  case  of  ne(o )  =  nm(o)  =  2  •  1012  cm-3  which  is  larger  than  the  steady 
state  electron  density  neo.  ne  then  grows  indefinitely  in  time,  verifying  the  instability.  One 
finds  that  Figs.  3.S  and  3.9  represent  the  general  time  evolution  of  the  discharge  in  the 
stiff  voltage  case,  depending  on  whether  ne(o )  is  larger  or  smaller  than  ne0. 

At  this  point  it  is  appropriate  to  proceed  to  the  numerical  simulation  of  the  finite 
external  impedance  discharge  case.  Again  setting  Va  =  9  kV/cm,  1  =  1  cm.  r  =  1  ^sec 
and  (3  =  10J  sec-1  =  Const(f).  Nevertheless,  in  contrast  to  the  stiff  voltage  case,  the  rate 
Eqs.  (3.5)  and  (3.6)  have  now  been  solved  simultaneously  with  the  circuit  Eq.  (3.14). 
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TIME  (nsec) 


Figure  3.8:  The  time  evolution  of  the  electron  and  metastable  densities  in  a  stiff  voltage 

discharge.  The  initial  electron  density  (n^O)  =  5xl011  cm*3)  is  lower  than  the  steady 
state  density  (neo  -  1.07xl012  cm-3)  at  the  given  voltage  (Vo  =  8.45  kV). 
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Figure  3.9:  The  time  evolution  of  the  electron  and  metastable  densities  in  a  stiff  voltage 

discharge  with  attachment  Initial  electron  density  (2xl012  cm'3)  is  higher  than  the 
steady  state  density  (n^  =  1.07xl012  cm-3)  at  the  given  voltage  (V0  =  8.45  kV). 
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The  code  also  makes  use  of  the  dependencies  of  the  rate  constants  shown  in  Figs.  3.4- 
3.6.  The  predictions  of  the  code  are  plotted  in  Figs.  3.10,  3.11  and  3.12.  Figure  3.10 
shows  the  temporal  variation  of  ne  and  E  for  the  circuit  impedance  of  RS  =  400  Qcm2. 
The  initial  electron  density  was  assumed  to  be  1012cm-3.  One  observes  that  the  electron 
density  increases  by  about  10%  to  1.07  •  1012cm-3  and  the  electric  field  decreases  from 
its  initial  value  of  9  kV/cm  to  S.45  kV/cm.  These  axe  the  steady  state  electron  density 
and  the  electric  field  under  the  above  conditions  (see  the  graphical  solution  in  Fig.  3.7). 
Figure  3.11  shows  the  effect  of  lowering  the  external  impedance  to  RS  =  10  cm2.  Under 
these  conditions,  as  can  be  seen  in  the  figure,  the  discharge  has  stiff  voltage  characteristics 
during  the  first  ~  10  //sec  of  the  discharge.  It  is  unstable  as  predicted  by  the  analysis  and 
the  electron  density  suddenly  increases  by  almost  three  orders  of  magnitude  after  10  /isec. 
The  discharge  becomes  stable  at  this  large  value  of  ne  ~  8.1014cm~3.  However,  in  practice, 
any  initial  spatial  nonuniformities  in  ne  amplified  in  time  by  the  instability  could  result 
in  streamer  formation  during  the  fast  transition  period.  We  shall  consider  the  details  of 
such  spatial  ir«*  ^Hties  in  the  next  chapter  of  this  report.  Finally,  Fig.  3.12  illustrates 
a  typically  situation  in  a  discharge  with  an  above-critical  external  impedance  (R  >  Rcr, 
see  the  discussion  in  the  previous  section).  In  the  case  considered  (see  Fig.  3.7)  Rcr  ~  520 
Qcm2,  while  the  results  presented  in  Fig.  3.12  correspond  to  the  case  R  =  S00  ficm2.  No 
steady  state  exists  in  this  case,  as  predicted  by  the  theory,  and  one  observes  in  Fig.  3.12 
that  the  discharge  is  quenched. 

3.5  CONCLUSIONS 

(i)  We  have  considered  the  volumetric  stability  dynamics  of  discharges  dominated  by 
two  step  Ionization.  In  the  presence  of  an  attacher  these  discharges  are  shown 
to  be  potentially  unstable  on  much  shorter  time  scales  than  those  characteristic  of 
conventional  thermal  instabilities. 

(ii)  The  theoretical  modeling  of  the  discharge  included  the  discharge  kinetics  coupled  to 
the  external  electric  circuit  and  its  response  to  varying  current  conditions.  A 
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Figure  3.1 1:  Curves  of  ne,  and  E  for  a  low  external  impedance  discharge  (RS  =  1  W  cm2). 

The  parameters  are  ne(0)  =10^  cn r3,  Vq  =  9  kV/cm,  L  =  1cm,  t  =  1msec  and  b  =  1CP 
sec-1.  The  discharge  is  unstable  during  the  inidal  evolution  phase  (t  <  15  msec). 

3-21 


SCIENCE  RESEARCH  LABORATORY 


Figure  3.12:  The  curves  of  ne  and  E  for  a  discharge  with  an  above-  critical  impedance 

(RS  =  800  Q  cm2).  All  other  parameters  are  the  same  as  in  Figs.  3.10  and  3.11.  No  steady 
state  exists  for  RS  >  RcS  =  520  Q  cm2  and  the  discharge  is  quenched. 
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number  of  discharge  configurations,  such  as  stiff  voltage,  stiff  current  and  finite  external 
impedance  were  considered  in  searching  for  a  steady  state  stable  regime  of  operation. 

(iii)  In  the  stiff  voltage  discharge  [the  electric  field  in  the  discharge  E  =  const  (t)]  the  steady 
state  exists  for  E  >  E*,  where  E*  is  the  value  at  which  the  characteristic  function  Q(E ) 
[see  Eq.  (3.9)]  vanishes.  Nevertheless,  even  if  the  steady  state  is  reached,  the  discharge 
is  volumetrically  unstable  under  the  stiff  voltage  condition.  Practically,  if  at  t  =  0 
the  initial  electron  density  ne(0)  >  nt0  ( ne0  being  the  steady  state  density)  the  current 
density  j  increases  indefinitely  on  a  finite  time  scale,  while  for  ne(0)  <  nt0  the  discharge 
is  quenched  (j(+)  —*  0).  The  linear  instability  temporal  rate  is  typically  proportional 
to  the  attachment  rate  0  (if  0  is  sufficiently  small). 

(iv)  The  stiff  current  discharge  [j  =  const(t)]  has  a  steady  state  at  all  currents  and  the 
discharge  is  volumetrically  stable,  provided  the  characteristic  function  F(E)  [see  Eq. 
(3.11)]  is  decreasing  function  of  the  electric  field. 

(v)  In  the  finite  external  impedance  discharge  case  we  found  the  existence  of  a  critical 
impedance  value  Rcr  such  that  there  exist  two  steady  state  when  0  <  R  <  Rcr  and  no 
steady  state  is  possible  when  R  >  Rcr.  Furthermore,  only  if  Rcr  is  large  enough,  the 
above  mentioned  steady  state  solutions  are  stable  in  a  part  of  the  interval  [0,  J?cr]. 

(vi)  The  theoretical  predictions  were  illustrated  in  the  example  of  an  atmospheric  C02  laser 
discharge  model.  The  numerical  solutions  of  the  kinetic  and  external  circuit  equations, 
underlying  the  problem  in  this  case  were  found  in  an  agreement  with  the  theory. 
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CHAPTER  4 


RF  STABILIZATION  OF  MOLECULAR  DISCHARGES 

4.1  INDUCTIVE  STABILIZATION  OF  THE  DISCHARGE  AGAINST 
STREAMER  FORMATION 

In  the  preceding  chapter,  the  volumetric  or  bulk  stability  of  the  discharge  has  been 
analyzed  and  stable  operating  regions  for  the  CO2  laser  discharge  have  been  identified. 
Although  a  discharge  can  be  volumetrically-stable,  micro-instabilities  can  lead  to  constric¬ 
tion  of  the  discharge  and  streamer  formation.  To  stabilize  the  discharge  against  streamer 
formation,  SRL  proposes  to  use  an  RF  pulsed  power  source.  Since  a  streamer  is  far  more 
inductive  than  the  discharge,  its  formation  and  growth  will  be  suppressed  because  the  RF 
power  component  will  be  excluded  from  the  region  of  a  localized  instability  as  it  attempts 
to  form.  The  inductance  of  the  localized  instability  region  is  responsible  for  excluding  RF 
currents  (and  therefore  power)  from  this  region. 

The  effectiveness  of  inductive  discharge  stabilization  can  be  evaluated  from  an  analysis 
of  the  equivalent  circuit  shown  in  Fig.  4.1.  The  discharge  impedance,  which  has  inductive 
and  resistive  components,  Lp  and  Rp,  while  operating  stably,  is  driven  by  an  oscillating 
current  source.  The  currents  through  these  components  are  Ii  and  In  respectively.  The 
voltage  drop  across  the  inductor  and  resistor  must  be  equal  hence 


iujLIi  =  R  In 


(4.1) 


By  conservation  of  current, 


h  + In  =  I 


(4.2) 


where  I  is  the  source  current.  Equations  (4.1)  and  (4.2)  can  be  solved  simultaneously  to 
give 


and 


R(R  —  iuL) 


(4.3) 
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of  discharge  by  inductance 
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For  efficient  discharge  pumping  the  inductive  current  should  be  minimized,  i.e.,  R  <<  uL, 
and  I  ~  In.  If  the  reverse  is  true,  then  I  ~  II  and  the  transfer  of  energy  to  the  discharge 
will  be  inefficient  especially  if  additional  resistive  losses  are  present  in  the  system. 

Should  the  discharge  begin  to  constrict,  the  inductance  will  increase  and  the  resistance 
will  decrease.  Under  these  conditions  the  frequency  of  RF  modulation  can  be  chosen 
to  exclude  RF  current  and  power  from  the  constricting  region  and  thereby  inhibit  arc 
formation. 

From  the  equations  above,  it  is  clear  that  the  RF  frequency,  f,  should  be  selected  such 
that  R£>/2tLd  >  f  >  Rs/2zLs  where  R,  and  L,  are  the  resistance  and  inductance  of 
the  perturbed  discharge  region.  These  constraints  on  the  modulation  frequency  have  a 
simple  physical  explanation.  The  constraint,  Rq/2-kLo  >  /,  insures  that  the  skin  depth 
at  the  modulation  frequency  is  larger  than  the  transverse  dimension  of  the  discharge.  This 
constraint  leads  to  spatially  uniform  discharge  pumping  of  the  laser  medium.  On  the 
other  hand,  the  constraint  R,/2ttL,  <  /  insures  that  the  skin  depth  at  the  modulation 
frequency  is  smaller  than  the  transverse  dimension  of  the  discharge  constriction.  Therefore, 
oscillating  discharge  power  is  prevented  from  penetrating  to  the  center  of  the  constricted 
region  and  continuing  the  process  of  streamer  formation.  Regions  of  ionization  instability 
in  CO2  laser  discharges  are  dominated  by  the  multistep  ionization  processes  discussed  in 
the  preceding  sec  Lion.  In  this  case,  the  ionization  rate,  t/,,  is  not  simply  a  function  of 
the  effective  electric  field,  E/N,  but  is  a  function  of  the  power  density  deposited  into  the 
discharge  E  •  J/N.  Consequently  the  limited  power  flow  into  the  constricted  discharge 
region  reduces  the  ionization  rate  in  the  region  and  enhances  stability. 

The  intrinsic  inductance  of  the  discharge  can  be  estimated  from  the  geometrical  con¬ 
siderations  shown  schematically  in  Fig.  4.2.  Also  shown  in  Fig.  4.2  is  the  shape  of  the 
magnetic  field  for  the  symmetric  current  condition  with  no  electron  density  perturbations. 
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The  discharge  inductance  is  defined  a s  Lp  =  j B  ■  ds/I.  For  a  square  aspect  ratio  (i.e.,  t 


=  h  in  Fig.  4.2) 


(4.5) 


where  L  is  the  length  of  the  active  volume  in  the  lasing  direction  and  fi0  is  the  permeability 
of  free  space. 

The  inductance  of  the  streamer  representing  the  electron  density  perturbation  can 
be  calculated  in  a  similar  maimer.  The  geometry  of  the  streamer  is  shown  in  Fig.  4.3. 
Assuming  that  the  streamer  has  collapsed  to  a  radius  a, 


,  not,  2 e 
L,  s  —In— 
2n  a 


(4.6) 


Dividing  Eq.  (4.6)  by  Eq.  (4.5)  one  gets  Ls/Lp  ~  2L(ln2£/a)/irL  Typically  £  =  0.15 
meters,  tn  2£/&  ~  4  and  L,/Lp  ~  20. 

To  calculate  the  required  frequency  of  modulation  the  resistance  of  the  CO2  discharge 
must  be  estimated.  The  conductivity  of  the  discharge  is  typically  3  x  10-2  mhos/m.  For 
a  discharge  cross  section  of  1.5  x  103  cm2  and  a  discharge  length,  £,  of  0.15  meters,  the 
discharge  resistance  is  35  ohms.  The  inductance  as  given  by  Eq.  (4.5)  is  7  nH.  So  for  the 
example  cited  above,  the  RF  modulation  frequency  should  be  in  the  range 

5.3  x  10 ‘^-Hz  <  /  <  8  x  10 sHz  (4.7) 

As  stated  above,  the  upper  limit  on  modulation  frequency,  800  MHz,  insures  that  the  skin 
depth  of  the  discharge  is  larger  than  the  transverse  dimension  of  the  discharge  so  that 
the  laser  medium  is  pumped  uniformly  and  efficiently.  The  remaining  inequality  relates 
to  inductive  stabilization  of  constrictions  in  the  discharge  current.  This  inequality  insures 
that  the  oscillating  portion  of  the  discharge  current  does  not  contribute  to  the  current  in 
the  constricted  region.  The  utility  of  this  inequality  is  best  understood  by  rewriting  it 
as  Rs/Rp  <  f(MHz)/ 53.  This  inequality  relates  the  streamer  resistance  (normalized  to 
the  stable  discharge  resistance)  to  the  modulation  frequency  required  to  insure  that  the 
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Figure  4.3:  Inductance  of  constricted  discharge 
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oscillating  current  does  not  flow  through  the  streamer.  Initially,  in  the  stable  discharge 
regime  Rs  is  infinite  and  therefore  RJRd  >>  // 53  MHz.  As  the  discharge  constricts, 
R,/Rd  falls  below  f/53  MHz  and  the  oscillating  current  is  cut  off  from  the  constricted 
region.  At  this  point,  the  excess  electron  density  in  the  constricted  region  should 'rapidly 
decay  thereby  driving  the  discharge  back  to  stable  operation.  If  the  modulation  frequency 
in  this  example  were  chosen  to  be  50  MHz.  then  stabilization  should  begin  to  occur  when 
the  streamer  resistance  begins  to  drop  below  the  discharge  resistance,  i.e.  as  the  streamer 
begins  to  form. 

4.2  TWO  DIMENSIONAL  RF  STABILITY  THEORY 

Two  rate  equations  are  important  in  determining  the  stability  of  the  discharge.  The 
first  is  the  electron  production  and  loss 

~  =  vnmne  -  an2t  -  (3ne  (4.S) 

at 

and  the  second  is  the  metastable  production  and  loss 

—p  =  (<7u)ncna  -  Tim/rm  (4.9) 

at 

where  ne  is  the  electron  density,  nm  is  the  metastable  density,  a  is  the  recombination  rate, 
/3  is  the  attachment  rate,  ( av )  is  the  electron  impact  metastable  production  rates  constant 
and  Tm  is  the  metastable  lifetime.  The  ionization  rate  constant  u  in  Eq.  (4.8)  is  assumed 
to  be  the  result  of  metastable  ionization.  Since  the  metastable  levels  have  a  much  smaller 
ionization  energy  than  the  ground  state  the  ionization  rate  of  the  discharge  is  dominated 
by  electron  impact  ionization  of  the  metastables.  The  stability  of  Eqs.  (4.S)  and  (4.9) 
will  be  analyzed  assuming  that  the  applied  discharge  is  sinusoidally  varying.  Hence  the 
various  rates  are  determined  by  the  local  RF  electric  field  of  frequency  w  which  is  generated 
externally  by  the  system  shown  in  Fig.  4.4.  From  the  analysis  performed  in  Phase  I  Science 
Research  Laboratory  (SRL)  has  shown  that  a  current  source  is  volumetrically  stable  and 
so  the  initial  analysis  will  assume  an  RF  current  source. 
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Figure  4.4:  Schematic  of  a  four  dimensional  discharge  and  electrical  circuit 
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We  shall  assume  now  that  the  rates  describing  the  variation  of  ne  and  nm  in  Eqs. 
(4.8)  and  (4.9)  (i.e.  vnm,  ane,  f3,  7 ne/nm)  are  much  smaller  than  w/27r,  so  that  one 
can  use  time  averaged  rate  constants  in  these  equations.  On  the  other  hand,  we  shall 
assume  that  u/2n  is  much  smaller  than  the  electron  energy  distribution  equilibration  rate 
(u ineiastic) ,  and  therefore  in  the  above  mentioned  time  averaging  process  one  can  use  the 
rates  corresponding  to  the  stationary  discharge  with  the  electric  field  equal  everywhere  to 
the  field  in  the  RF  discharge  at  a  given  time.  Thus  we  assume 


unm  «  ui/2n  «  vinela3tic  (4.10) 

Assuming  a  100  Torr  discharge,  vnm  ~  106  sec-1  and  i/‘*w  =  10u  sec-1,  so  that  1  MHz 
<<  /  =  w/27r  <<  102  GHz. 

At  this  stage  it  will  also  be  assumed  that  the  vacuum  wavelength  of  our  oscillator 
(A  =  2irc/u; )  is  much  larger  than  L,  so  that  in  the  vicinity  of  the  discharge  one  can 
neglect  the  relativistic  retardation  of  time  due  to  the  extended  volume,  and  treat  the 
electromagnetic  problem  as  being  quasistationary.  For  i  =  lm,  the  condition  L  «  A 
yields  /  <<  c/L  =  3  •  108  sec-1.  So  that  we  are  limiting  our  discussion  to  frequencies  in 
the  range  10  MHz  <  /  <  100  MHz. 

The  quasistationarity  condition  means  that  in  the  Maxwell  equation  for  Vx  S  one  can 
neglect  the  displacement  current  and  use  the  “magnetstatic”  law 

Vx  9  =  t  ^ (t)  (4-n) 

while  E  is  given  by  the  Faraday  law  and  the  quasi-  neutrality 


Vx  E= 


Id  B 
c  dt 


(4.12) 


V-  E=  0 


Equations  (4.11),  (4.12)  and  (4.13)  yield 


V2  E=  9  J 


dt 
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where  the  current  density  in  the  plasma  is  given  by  j  =  a  E ,  with  the  conductivity  a  — 


—ent/j.  The  simplest  model  for  the  mobility  /x  will  be  used  to  give 


mvcou 


(4.15) 


Now  we  shall  consider  the  steady  state  of  Eqs.  (4.8)  -  (4.9)  assuming  cylindrical 


symmetry  and  that 


E=  Re(E0e~ut)et 


Therefore,  all  the  parameters  of  the  steady  state  (neo,  nmo,  <r0,E0 )  are  functions  of  r  alone. 


In  steady  state  one  obtains 


—  yTTleo 


wyr  —  a 


All  the  rate  constants  are  temporally  averaged  over  the  voltage  oscillation. 

Returning  back  to  Eq.  (4.7),  the  steady  state  equation  gives 

1  d  (  dE„\  Ane2neo  uj  „ 

(  r—  +  i - - - E0  =  0 

r  dr  \  dr  J  me1  ucou 

The  solution  of  this  equation  is  (assuming  neo  ~  const(r)) 


(4.16) 

(4.1?) 


(4.  IS) 


=  AJJxr) 


(4.10) 


where 


2  _  .4?r e2nt0  u  _  2 i 


X  =  * 


and  S  is  the  collisional  skin-depth 


me2  vcoll  62 


_  /  emucou  \ 

\2Ke2neou)J 


(4.20) 


(4.21) 


The  constant  .4  in  Eq.  (4.19)  can  be  found  from  the  condition  that  the  magnetic  field  on 
the  surface  of  the  plasma  equals  to  2 10/ciZ,  where  I0  =  V0/Z  is  the  current  given  by  the 


external  current  source.  From  the  Faraday  law 


^  =  (Vxi),  =  -^ 


(4.22) 
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(4.23) 


- - xAJiixr) 

c J 


and  therfore 


A  =  i 


2 10UJ 


c2Rxj 


In  the  limit  of  low  frequencies,  when 


R  ,  c2  mu  con 

7  «  1  W  <<C  2Tre2nt0R2 


(4.24) 


(4.25) 


we  can  use  the  small  argument  expansion  of  the  Bessel  function  to  obtain  the  expression 
for  E0 : 


2I,u>  (  S  V  r  i  fr\ 

E°{r)-  —  \J)  I'-zU) 


(4.26) 


Condition  (4.25)  is  important,  since,  as  follows  from  (4.26),  it  allows  one  to  have  a  radially 
uniform  discharge  in  the  steady  state.  For  example,  in  the  case  of  a  100  Torr  laser  discharge 
{vcoii  ~  3-  1011  sec-1),  R  =  3  cm  and  neo  =  1012  cm-3,  (4.25)  yields  u  «  2-1010  (/  <<  3 
GHz)  which  is  within  the  quasistationary  limit  (/  <  100  MHz)  described  earlier. 

When  (4.25)  is  satisfied,  Eq.  (4.26)  yields 


Eo(r)  ~ 


2 10u 


f£V~-= 

\RJ  nt0 


const(r) 


It  is  precisely  this  inverse  dependence  of  E0  on  ne  which  assures  the  global  stability  of 
the  discharge  in  the  current  source  scheme.  Any  increase  (or  decrease)  of  the  total  density 
leads  to  the  decrease  (or  increase)  of  E0  and  therfore  to  the  decrease  (or  increase)  of  the 
electron  and  metastable  production  rates.  This  suppresses  the  initial  perturbation.  The 
details  of  such  a  global  stabilization  by  the  current  source  are  similar  to  those  already 
considered  previously  for  the  DC  discharge.  In  contrast  to  the  DC  discharge,  however, 
the  RF  current  source  also  allows  to  control  the  local,  streamer  type  instabilities  via  the 
skin-effect.  This  analysis  will  be  presented  next. 


4-11 


SCIENCE  RESEARCH  LABORATORY 


4.3  STABILITY  ANALYSIS 


Assume  that  the  large  skin-depth  condition  (4.25)  is  satisfied  for  a  given  steady  state 
so  that  ne0,  rcmo,  <r0,  E0  are  "almost”  (to  order  (R/6)2)  uniform  throughout  the  plasma 
volume.  One  can  introduce  a  small  local  streamer  type  density  perturbation  of  a  char¬ 
acteristic  transverse  size  A.  The  diffusion  term  will  be  added  to  the  rate  equations  since 
non- uniformities  are  of  interest.  Hence  one  can  write 

=jnmnt-  an]-  /3ne  +  DaV2ne  (4.27) 

Ot 

—  7ne  —  +  Dmv2nm  (4.2S) 

ot  T 

where  Da  and  Dm  are  the  ambipolar  and  metastable  atomic  diffusion  coefficients  respec¬ 
tively.  Dm  is  of  the  order  l/2vth.\  ~  102  cm2/sec/Torr  while  Da  may  probably  be  10  times 
larger  [Da  ~  103  cm2/sec  Torr). 

The  perturabation  in  ne  leads  to  a  local  conductivity  perturbation  which  in  turn  results 
in  the  decrease  or  increase  in  the  (rms)  electric  field  value,  depending  on  the  sign  of  the 
density  perturbation.  The  effect  is  similar  to  the  skin-effect  and  can  be  used  to  suppress 
a  possible  absolute  instability  in  the  discharge  (as  that  associated  with  the  attachment). 
This  stabilization  phenomenon  will  be  considered  in  more  details.  In  what  follows,  the 
diffusion  of  metastables  will  be  neglected  for  simplicity.  It  will  also  be  assumed  that  r_1 
is  the  fastest  rate  in  Eqs.  (4.27)  and  (4.28),  i.e., 

-  >>  immo,  an„,  /?,  Da/ A2  (4-29) 

T 

Typically  the  rate  r_1  is  in  the  few  MHz  region,  while  the  other  rates  belong  to  the 
sub-MHz  scale.  For  times  t  >  l//s,  condition  (4.29)  Eq.  (4.28)  to  be  replaced  by 

nm(t)  ~  r~(ne(t)  (4.30) 


(this  is  the  fast  scale  analysis  of  the  previous  study)  Then  Eq.  (4.27)  becomes 


=  (t'7r  ~  » )n\  -  f3nt  +  DaV2ne 
at 
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If  ne  =  const,  Eq.  (4.31)  again  yields  the  steady  state  Eq.  (4.17). 

At  this  point  the  perturbation  analysis,  will  be  employed  i.e.  solution  of  Eqs.  (4.14) 


and  (4.31)  in  the  form 


Mg  —  4"  ft*  1 


E  =  (E0  +  E1)e,uteJ 

where  all  the  perturbations  (nel,£ i)  vary  slowly  in  time 


(4.32) 


1  dln£ 

tjj  dt 


«  l,f  =  nel,  E] 


The  perturbed  equations  become  (after  the  linearization): 


2  t-i  •  47re2cj  _  47re2  / 9nei  d£i 

V2£x  =  - (nel£0  +  neoEx)  +  — - (  £0  +  nt0— — 

c^mv-oii  czmuC0u  \  at  at 


(4.33) 


=  [2(1/77-  -  a)neo  -  (3 ]  n*i  +  DaV2nei  +  [(i/7r  -  a)'n20  -  /?'ne0]  £1  (4.34) 

Here  (...)'  =  d/dEa  (...).  Finally,  neglecting  the  selected  term  in  Eq.  (4.33)  and  using  the 
steady  state  relation  (4.17)  in  rewriting  Eq.  (4.34)  one  obtains 


(4.35) 


(4.36) 


=  &nti  +  £aV2nei  +  /?neo  _  £]  £j  (4.36) 

dt  [{u~fT  —  a)  p. 

The  skin-depth  8  has  again  been  introduced  in  Eq.  (4.35). 

At  this  stage  a  solution  of  Eqs.  (4.35)-(4.36)  in  the  form  ntl,Ei  ~  e,k  r  are  searched 
for,  where  |£|2  is  of  0(2~/A)  Then  one  has 


1  2  1  “  r>  “  J~‘° 

-k  +  —  E 1  = -— - nt 


where 


=  (3  ~  Dak2)nel  +  '3neoB'E\ 


B'  -  ^7T  Z  Q Y  _  £. 

{v~fT  -  a)  3 


By  substituting  E\  from  (4.37)  into  (4.3S)  one  obtains 

dneX  uo  n  ,2  4 3neoEnB\  ,  . 23neoEnB'k26 2, 

—  =  ^-D‘k  -  W  +  4  )  +  ,_  W+4  -jn'' 


(4.3S) 


(4.39) 


(4.40) 
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Since  k  >  2t/R  and  S  >>  R  one  finds  (fca)4  >>  4  and  therefore  (4.40)  reduces  to 


d(lnnei)  2  4/3nco£’0J3'  .  r  2 (3ne0E0B' 


=  f3~  Dak*  - 


+  L 


dt  r  fc262 

The  last  equation  shows  that  a  stable  operation  corresponds  to  the  case 


(4.41) 


$  =  /3  -  Dak2  - 


40ncoEoB' 


<  0 


(4.42) 


One  observes  that  there  exist  two  stabilizing  phenomena  in  the  system.  The  first  is 
due  to  the  diffusion  and  is  effective  in  stabilizing  short  scale  density  perturbations  The 
diffusion  simply  smoothes  the  short  scale  non- uniformities.  The  second  phenomenum  is 
basically  the  skin  effect.  It  is  more  important  in  suppressing  larger  scale  perturbations, 
espesially  those  of  size  comparable  to  the  radious  of  the  plasma  column. 

The  interplay  between  the  two  stabilizing  processes  can  be  demonstrated  by  drawing 
the  curve  showing  the  dependence  of  $  =  $(fc2)  (see  Eq.  (4.42)).  A  typical  situation  is 
shown  in  Fig.  4.5. 

${k2)  reaches  its  maximum  value 


0-2.05 


1/3 


at 


2  _  /8 Pnt0E0B' 
m  \  Da6* 


1/3 


(4.43) 


(4.44) 


Therefore,  in  case  km  >  2ir / R.  the  stability  condition  for  the  entire  k  -  space  becomes 

(4.45) 


sG(_n^ 


Note  that  the  stable  operation  is  easier  at  higher  frequencies  and  plasma  densities. 
Note  also  that  the  final  value  of  Da  is  important,  since  (38)  can  not  be  satisified  for 

Da—+  0. 
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Figure  4.5:  Plot  showing  the  functional  dependence  of  <t»  (k2)  vs  k2 
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4.4  NUMERICAL  SIMULATIONS  OF  THE  TWO  DIMENSIONAL  CODE 

SRL  has  developed  a  code  based  on  the  two  dimensional  theory  presented  earlier  in 
this  section.  The  basic  equation  is  (4.18)  which  can  be  in  rectlinear  coordinates  as 

d2E  i47re2n,„  u 


+ 


dx 2  '  me 2  ucou  0 


Ea  =  0 


(4.46) 


subject  to  the  boundary  conditions 


dE, 


dE. 

dx  'I=0  ”  °  and  dx  ll=t 


2ujL 


where  I„  is  the  total  current  per  unit  length.  Equation  (4.46)  is  valid  for  a  slab  geometry 
i.e.  where  the  z  dimension  is  much  longer  them  x-dimension.  Such  a  geometry  conforms 
more  to  the  CO2  lasers  than  the  radial  geometry.  As  discussed  in  Chapter  2  the  frequency 
of  osicllation  /  >>  ane  hence  the  electron  density  will  remain  essentially  constant  during 
a  period  of  oscillation.  This  allows  one  to  use  time  averaged  rates  for  o,  u  and  ( av ).  These 
time  averaged  rates  are  defined  by  the  following  integral 

t/uj  r 

£  =  —  j  £(|  E  |  sinu >t)dt  =  —  J £(E0sins)ds  =  f(|  E  |)  (4-47) 

O  O 

where  £  =  1/,  a  and  {crv). 

In  Eq.  (4.47)  |  E  j  is  just  the  amplitude  of  the  electric  field 


E  |=  {(M)2  +  M2}1/2 


(4.48) 


For  efficient  simulation  one  has  to  obtain  analytic  expressions  for  P,  a  and  (av)  as  func¬ 
tions  of  j  E  |.  It  is  also  convenient  to  introduce  a  dimensionless  coordiante  p  =  x/w  in 
which  case  Eq.  (4.46)  can  be  rewritten  as 


d“  /w\‘i 

opE  +  2'(j)  £  =  ° 


(4.49) 


where  the  skin  depth  6  is  given  by 

8 = (  c  y/2  =  ( cmuc°i1  \ 

\2irauj)  \27re2nea;  / 


1/2 
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As  stated  in  Chapter  2  the  discharge  width  u>  <<  6  and  so  one  can  solve  Eq_  (4.46) 
by  a  perturbation  expansion  by  introducing  the  small  parameter  A  =  (w/8)2 .  The  electric 
field  can  be  expanded  as 

E  =  E0  +  E\  4-  Ei . 

where  the  terms  are  ordered  in  powers  of  A  then  Eq.  (4.46)  gives  the  following  orders 


_  0 
dp 2 

(4.50) 

L  +  2iAE0  =  0 

(4.51) 

2 

L  +  2iAEi  =  0 

(4.52) 

Equation  (4.50)  perdicts  Eq  is  constant  which  after  substitution  into  (4.51)  results  in 

p  p' 

Ei(p,t)  =  —2iE0  J  dp'  J  A {p")dp"  (4.53) 


Similarly 


E,{p,t)  =  -2 i  j  dp'  J'  dpE^&dp" 

O 

etc.  The  SRL  code  solves  the  first  order  equation  given  by 

E(p,t)  =  E0(l-  2  i  j'  dp'  j  A  (/,  t)dp" 


(4.54) 


(4.55) 


The  constant  E0  can  be  found  from  the  boundary  condition  on  the  total  current 

l 


dp  l»-' 


=  -2 iE.J  &(p',t)p'dp' 


E.=  ^./  J\(p',t)p'dp'  (4.50) 

The  numerical  code  results  will  be  presented  next.  The  particular  slab  geometry  analysed 
is  shown  in  Fig.  4.6.  The  discharge  is  assumed  to  be  2.4  cm  wide  in  the  x-  direction  and 
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X 


Schematic  of  discharge  geometry 
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long  in  the  z  direction.  The  electric  field  is  applied  in  the  y  direction.  The  origin  of  the 
abscissa  is  located  at  the  center  of  the  discharge  i.e.  the  discharge  extends  between  -1.2 
cm  to  1.2  cm  in  the  x-  direction.  The  stiff  current  drive  will  be  discussed  first  since  it  was 
the  most  stable. 

4.4.1  A  Stiff  Current  Case 

Figure  4.7  shows  the  temporal  and  spatial  variation  of  nt  for  the  discharge  between 
0  <  x  <  1.2  cm.  The  discharge  is  assumed  to  be  symmetric  hence  —  1.2  cm  <  x  <  0  will 
just  be  the  mirror  image  about  the  ordinate  of  the  plots  shown  in  Fig.  4.7.  Also  it  was 
assumed  that  the  total  current  was  2.4  A/cm  and  the  applied  electric  field  was  400  V/cm 
and  the  CO2  laser  mixture  was  3/2/1  at  a  total  pressure  of  30  Torr.  The  initial  spatial 
electron  density  was  assumed  to  be  slightly  greater  than  1011  cm-3  with  a  non  uniformity 
of  5%  between  the  centre  of  the  discharge  and  the  edge.  The  attachment  rate  for  electrons 
was  assumed  to  be  50  /is.  The  spatial  variation  of  the  electrons  is  shown  every  5  /is.  The 
initial  spatial  profile  is  labled  1.  The  curve  labled  2  is  spatial  profile  5  /is  later  etc.  The 
electron  density  increased  by  20%  in  the  first  10  /is  and  then  decrease  by  10%  during  the 
next  50  /is.  Also  the  non  uniformity  at  the  end  of  45  /is  increased  to  about  10%.  Figure 
4.S  shows  the  results  of  the  next  45  /is.  The  curves  in  Figure  4.8  are  spaced  every  10 
/isec.  The  curve  labled  1  in  Fig  4.8  is  the  same  as  the  curve  labled  10  in  Fig.  4.7.  The 
curve  labled  2  shows  the  spatial  variation  of  ne  10/is  later,  etc.  The  dominant  effect  is  the 
continual  increase  in  the  non  uniformity.  The  uniformity  at  the  end  of  85  /is  is  25%. 

The  result  shown  in  Fig.  4.7  and  4.8  show  for  the  conditions  assumed  the  discharge 
was  stable  for  about  100  /is.  However  the  spatial  non- uniformity  increases  temporally 
indicating  the  discharge  will  eventually  be  unstable  and  arcing  will  occur.  From  the  many 
code  runs  analysed  it  is  clear  that  a  spatial  non-uniformity  and  a  finite  attachment  rate 
results  in  an  unstable  discharge  even  for  a  stiff  current  source  and  the  key  for  stable 
operation  is  to  obtain  as  uniform  an  initial  electron  density  as  possible  and  to  keep  the 
electron  attachment  rate  as  small  as  possible. 
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DISTANCE  FROM  CENTER  (cm) 


Figure  4.7:  Plot  showing  the  spatial  and  temporal  variation 
of  the  electron  density  ne  .  Curve  1  is  the  initial 
spatial  variation  of  ne  at  t  =  0.  Curve  2  is  the 
spatial  variation  of  ne  at  t  =  5  p.s,  etc.  The  power 
supply  is  assumed  to  be  a  stiff  current  source. 
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ELECTRON  DENSITY  AT  DIFFERENT  TIMES 


k  DISTANCE  FROM  CENTER 


4.4.2  Stiff  Voltage  Source 

Figures  4.9  and  4.10  show  the  results  of  a  stiff  voltage  source  for  the  same  conditions 
discussed  previously.  In  Fig.  4.9  the  initial  electron  density  was  again  chosen  to  be  about 
l.lxlO11  cm-3.  The  electron  density  increases  by  over  an  order  of  magnitude  in  30  pts. 
In  this  figure  the  spatial  variation  in  the  electron  density  is  shown  every  2  fis  as  apposed 
to  every  5  \i s  in  Fig.  4.7  and  every  10  /xs  in  Fig.  4.8.  This  discharge  is  clearly  unstable 
and  it  will  result  in  an  arc  well  before  the  100  \is  pulse  length.  Figure  4.10  shows  the 
results  of  initializing  the  electron  density  to  1010  cm-3.  In  this  case  the  discharge  is  clearly 
quenched.  These  results  are  remincent  of  the  one  dimensional  stiff  voltage  results  shown  in 
Figs.  3.11  and  3.12  where  discharge  is  unstable  and  it  proceeds  to  an  arc  or  its  is  quenched 
depending  on  the  initial  electron  density.  So  once  again  it  is  clear  that  the  stiff  current 
source  is  more  stable  than  a  stiff  voltage  source. 
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ELECTRON  DENSITY  AT  DIFFERENT  TIMES 


TIME 

TEMPORAL  SEPARATION 
BETWEEN  CURVES 
=  3  usees 


Figure  4.9:  Same  as  Figures  4.7  and  4.8  except  that  power  supply 
is  assumed  to  be  a  stiff  voltage  source  and  the  spadal 
variation  of  ne  is  shown  every  3  fis. 
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CHAPTER  5 


CONCEPTUAL  DESIGN  OF  A  REPETITIVELY  PULSED  C02  LASER 

In  this  chapter  a  conceptual  design  of  the  laser  head  and  the  RF  pulse  power  to  drive 
it  will  be  presented.  An  assembly  drawing  of  the  laser  head  is  shown  in  Fig.  5.1.  The 
laser  head  is  designed  with  several  spark  plugs  that  provide  UV  preionization.  The  waste 
heat  will  be  removed  by  conduction  to  the  side  walls.  Such  a  design  eliminates  the  flow 
loop  and  results  is  a  considerably  lighter  weight  C02  laser.  Also  the  design  is  simpler  and 
subsystems  such  as  a  blower  and  sidewall  mufflers  are  eliminated.  The  anode  and  cathode 
were  designed  using  a  Chang  profile  thereby  insuring  that  the  electric  field  at  the  edge  of 
the  electrodes  do  not  increase  substantially.  Such  an  increase  in  the  electric  field  could 
result  in  arcing.  The  motivation  behind  choosing  the  laser  head  design  will  be  discussed 
in  the  remainder  of  this  section. 

5.1  GEOMETRIC  CONFIGURATION  OF  THE  C02  LASER  HEAD 

In  order  to  size  and  configure  the  laser  head,  several  issues  must  be  considered  including 

•  Specific  extractable  laser  energy 

•  Gas  Cooling  Method 

•  Preionization  Method 

•  Efficiency  of  Optical  Extraction 

The  specific  energy  extractable  from  C02  lasers  is  typically  50  J /liter  atm.  SRL 
proposes  to  operate  the  laser  at  30  Torr  or  1/25  of  an  atmosphere  and  hence  we  plan 
to  extract  approximately  2  J /liter.  Therefore  a  6  liter  laser  volume  will  be  adequate  for 
demonstrating  the  12  J  laser  pulse  energy  required  for  SDIO  applications. 

The  next  design  consideration  centers  on  the  fact  that  one  of  the  dimensions  is  con¬ 
strained  by  conduction  cooling  to  be  about  2.4  cm.  This  constraint  comes  from  the  fact 
that  the  C02  laser  efficiency  decreases  at  elevated  temperatures  hence  there  is  a  maxi¬ 
mum  temperature  difference  that  can  be  tolerated  between  the  walls  and  the  axis  of  the 
discharge  cell  under  repetitive  pulse  operation. 
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For  efficient  laser  power  extraction  the  integrated  optical  length  through  the  laser 
device  must  be  >  500  cm.  SRL  proposes  to  use  a  cylindrically  expanding  beam  that  will 
be  folded  into  four  passes.  Therefore  a  discharge  length  of  125  cm  will  be  adequate.  The 
third  dimension  then  must  to  be  20  cm  to  provide  the  required  6  liter  volume.  Figure  5.2a 
schematically  shows  the  discharge  volume  and  Fig.  5.2b  shows  the  four  folds  that  enables 
one  to  take  the  expanding  beam  geometry  shown  in  Fig.  5.3  and  make  it  compatible  with 
a  rectangular  geometry.  The  discharge  voltage  will  be  applied  across  the  20  cm  length  and 
the  UV  preionization  will  be  introduced  by  the  use  of  spark  plugs  that  are  placed  behind 
a  cathode  mesh. 

5.2  C02  LASER  DISCHARGE  CHARACTERISTICS 

The  performance  of  C02  laser  discharges  is  determined  by  several  parameters  including 
electron  and  gas  temperature,  gas  mixture  and  the  type  and  level  of  preionization.  For 
efficient  laser  operation  one  must  choose  the  appropriate  laser  mixture  and  electric  field. 
The  mixture  of  choice  for  efficient  long  pulse  CO2  laser  operation  is  a  3/2/1,  He/N2/C02 
mixture.  The  rate  constants  for  the  ionization  and  excitation  of  the  electronic  states  of 
N2,  the  fraction  of  discharge  energy  that  goes  into  vibrational  excitation,  the  drift  velocity 
and  electron  temperature  as  predicted  by  the  Boltzmann  code  are  shown  in  Figs.  3.4. 
3.5  and  3.6.  From  Fig.  3.4,  in  which  Te.  vd  and  the  efficiency  of  exciting  the  vibrational 
levels  of  N2  are  plotted  against  the  discharge  electric  field,  it  is  clear  that  efficient  C02 
laser  operation  occurs  at  electric  fields  between  5-10  kV/cm-atm.  Thus  for  a  total  mixture 
pressure  of  30  Torr  the  appropriate  electric  field  will  be  200-400  V/cm.  A  consistent  set  of 
C'Oj  laser  discharge  characteristics  for  generating  a  12  joule  laser  pulse  (30  ^sec)  are  given 
in  Table  1.  These  discharge  characteristics  have  been  developed  using  the  comprehensive 
SRL  pulsed  C02  laser  model. 
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Figure  5.2a:  C02  gain  medium  geometry 
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Figure  5.3:  Expanding  beam  amplifier  concept 
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Table  1:  CO2  Laser  Discharge  Characteristics 


-  Discharge  Volume 

-  Discharge  Electric  Field 

-  Discharge  Voltage 

-  Current  Density 

-  Total  Current 

-  Discharge  Impedance 

-  Energy  Delivered  to  Discharge 


2.4  x  20  x  125  cm3 

400  V/cm 

8  kV 

1  A /cm2 

240  A 

33.3  ft 

72  J 


The  CO2  laser  to  be  designed  in  the  proposed  effort  will  be  UV  preionized.  Typically 
UV  preionization  results  in  a  low  level  of  electron  density  (~  106  -  108  cm-3)  that  has  to  be 
avalanche  four  to  six  orders  of  magnitude.  This  laser  will  also  have  a  prepulse  at  the  front 
end  of  the  main  discharge  pulse.  It  is  convenient  to  design  the  prepulse  to  have  a  voltage 
which  is  approximately  twice  the  voltage  of  the  main  discharge  pulse  which  means  that  the 
electric  field  will  be  about  600-S00  kV/cm  during  this  250  nsec  prepulse  duration.  Figure 
5.4  shows  the  variation  of  the  ground  state  ionization  rate  constant  va  as  a  function  of  the 
applied  electric  field.  At  an  electric  field  of  600  V/cm,  the  value  of  ua  is  about  2  x  10~10 
cm3/sec,  so  for  a  density  of  1018  cm-3,  the  ionization  rate  uana  ~  108  electrons/cm3 /sec 
and  the  electron  density  will  exponentiate  every  20  nsec.  As  discussed  earlier  in  this  section, 
the  electron  density  must  increase  4-6  orders  of  magnitude  or  9-14  exponentiations  which 
corresponds  to  a  prepulse  duration  of  200-300  nsec.  At  the  larger  electron  densities  of  1012 
cm-3  that  are  required  for  efficiently  pumping  CO2  lasers,  the  dominant  ionization  process 
will  be  two  step  ionization  from  metastable  levels.  In  this  ionization  process,  discharge 
electrons  first  produce  metastable  electronic  states  and  then  in  a  second  collision  take 
these  excited  states  into  the  ionization  continuum.  Under  SDIO/ONR  funding,  SRL  has 
developed  a  pulsed  CO2  laser  driver  which  efficiently  uses  a  single  circuit  to  provide  the 
preionizer  pulse,  the  prepulse,  and  a  main  discharge  pulse. 

5.3  OPTICAL  CONFIGURATION 

An  amplifier  is  much  more  desirable  than  an  oscillator  if  a  C02  laser  beam  optical 
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quality  of  better  than  1.3  XDL  and  a  spectral  width  of  less  than  100  kHz  are  required. 
However,  to  reduce  the  weight  of  the  system  the  amplifier  gain  must  be  as  large  as  pos¬ 
sible.  This  last  criterion  of  large  amplifier  gain  can  be  met  by  using  the  expanding  beam 
concept. W  This  design  concept  which  has  been  pioneered  by  and  experimentally  verified 
by  SRL  under  SDIO  sponsorship  makes  use  of  optical  beam  whose  area  changes  as  it  prop¬ 
agates  through  the  amplifier.  Such  a  configuration  is  shown  in  Fig.  5.3.  The  injected 
beam  which  has  a  cross-sectional  area  of  1  x  2.4  cm2  expands  cylindrical!}1-  to  an  area  of  9 
x  2.4  cm2  over  the  gain  length  of  500  cm.  From  this  geometry  it  is  clear  that  the  amplifier 
gain  is  the  product  of  the  area  gain  gA  and  active  media  gain  exp(gL),  where  gL  just  the 
laser  gain  length  product  along  a  ray  through  the  gain  medium.  Amplified  spontaneous 
emission  (ASE)  and  parasitic  modes  that  can  decrease  the  laser  efficiency  depend  on  (exp 
gL)  only.  Hence  this  geometry  lends  itself  to  a  much  higher  gain  per  stage  than  the  usual 
rectangular  extraction  geometry. 

This  fact  is  shown  graphically  in  Fig.  5.4  which  is  a  plot  of  the  amplifier  gain  vs 
extraction  efficiency  for  two  cases.  The  first  curve  is  for  the  expanding  beam  geometry 
shown  in  Fig.  5.4.  Note  as  the  amplifier  gain  decreases  (this  corresponds  to  increasing  the 
injection  level)  the  extracted  efficiency  increases.  Also  shown  in  Fig.  5.4  is  the  amplifier 
gain  vs  extraction  efficiency  for  the  rectangular  cavity  shown  in  Fig.  5.5  which  has  the 
volume  and  gain  length  as  the  expanding  beam  laser  cavity  of  Fig.  5.3.  The  extraction 
efficiency  of  the  rectangular  cavity  is  always  lower  than  the  expanding  beam.  Even  more 
significant,  however,  is  that  the  conventional  laser  cavity  design  will  need  isolation  for  a 
gain  of  greater  than  100.  This  becomes  clear  upon  looking  at  curve  4  that  predicts  the 
extraction  efficiency  of  a  rectangular  geometry  in  the  presence  of  1%  specular  reflection 
from  windows  or  other  optical  elements  in  the  beam  train.  The  rapid  decrease  of  the 
extraction  efficiency  for  stage  gains  of  >  100  is  a  result  of  amplified  spontaneous  emission 
(ASE).  From  curve  2  it  is  clear  that  the  EBL  is  relatively  insensitive  to  ASE  which  is 
due  to  the  beam  expansion  which  reduces  the  intensity  feedback  by  specular  reflection 
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Figure  5.4:  Curves  showing  the  extraction  efficiency  as  a  function  of  stage  gain  for  the 
EBL  and  the  rectangular  amplifiers. 
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by  an  order  of  magnitude.  These  ASE  calculations  axe  discussed  in  greater  detail  in  the 
following  subsection.  From  the  curves  in  Fig.  5.4  it  is  clear  that  one  cannot  build  an 
COj  laser  amplifier  as  a  single  500  cm  gain  length  for  gains  >  100  because  of  ASE  and 
parasitic  modes.  This  can  be  compared  to  an  amplifier  gain  of  1000  for  the  expanding 
beam  geometry. 

The  active  medium  as  shown  in  Fig.  5.6  is  unconventional  and  it  is  difficult  to  configure 
a  CO2  laser  discharge  into  such  a  volume.  However,  the  optical  beam  can  easily  be  folded 
and  into  a  rectangular  laser  discharge  volume  so  the  single  pass  configuration  shown  in 
Fig.  5.3  can  be  reduced  to  the  four  pass  rectangular  discharge  region  shown  in  Figs.  5.2b 
and  5.7.  The  length  of  the  discharge  region  is  reduced  to  125  cm.  The  cross-sectional  area 
is  20  x  2.4  cm'. 

5.4  LASER  GAIN  AND  EXTRACTION  EFFICIENCY 

Having  established  the  discharge  conditions  and  the  optical  configuration,  the  laser 
gain  and  optical  extraction  efficiency  can  be  calculated.  To  accomplish  this  calculation, 
SRL  has  used  a  combination  of  C02  laser  codes  which  have  been  developed  over  the  past 
4  years.  The  first  is  a  C02  laser/discharge  code  that  predicts  the  small  signal  gain  for  the 
discharge  conditions  given  in  Table  1.  The  small  signal  gain  as  predicted  by  that  code  is 
shown  by  the  curve  in  Fig.  5.7.  It  is  interesting  to  note  that  the  gain  increases  linearly  for 
the  first  10  ps  indicative  of  the  long  life  time  ~  100  p s  of  the  upper  level.  The  predicted 
laser  characteristics  are  summarized  in  Table  2. 


Table  2: 

! 


Laser  Characteristics 


•  Peak  small  signal  gain 

•  Injected  power 

•  Peak  output  power 

•  Energy  per  pulse 

• 


3.65  x  10  2  cm  1 
1  kW 
800  kW 
12  joules 
72  joules 


Discharge  energy  deposited 
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Figure  5.6:  Schematic  of  C02  laser  head 
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SMALL  SIGNAL  GAIN  (%  /  cm)  OUTPUT  LASER  POWER  (kw) 


TEMPORAL  PULSE  SHAPE 


•  Intrinsic  Efficiency 


16.7% 


To  calculate  the  optical  extraction  efficiency,  SRL  has  used  an  existing  expanding  beam 
code.  The  inputs  to  this  code  are  the  small  signal  gain  and  saturation  flux  as  predicted  by 
the  SRL  CO2  kinetics/laser  code  and  the  active  medium  geometry  shown  in  Fig.  5.7.  The 
output  laser  flux  as  predicted  by  the  code  is  shown  in  the  lower  curve  in  Fig.  5.7.  The 
peak  output  power  is  almost  S00  kW,  the  integrated  single  pulse  energy  is  12  joules  and 
the  overall  laser  power  extraction  efficiency  is  16.7%. 

5.5  AMPLIFIED  SPONTANEOUS  EMISSION  (ASE) 

The  stage  gain  of  a  CO2  laser  amplifier  will  be  determined  by  constraints  imposed  by 
ASE.  Hence  it  is  important  to  model  the  ASE  and  to  estimate  the  loss  of  efficiency  imposed 
by  this  process.  One  can  start  by  writing  the  intensity  radiated  from  an  elemental  volume 
in  the  laser 

A<£  =  N*/t  =  g<t>3  =  g0<t>s/(l  +  <t>/<t>s)  (5.1) 

ASE  is  just  the  amplification  of  A 6  along  a  given  direction.  The  longest  dimension,  which 
is  usually  the  laser  length  along  the  optical  axis,  will  result  in  the  largest  ASE  flux 

*ASE=i£~l  dr'urm7)exp{5(r-r')}  (5-2) 

In  Equation  (5.2)  Q  is  the  solid  angle  subtended  by  the  radiation.  From  Eq!  (5.2)  it  is 
clear  that  4>ase  increases  exponentially  as  the  gain-length  product  gL  increases.  Numerical 
integration  of  Eq.  (5.2)  indicates  that  efficient  laser  power  extraction  from  conventional 
CO2  laser  single  pass  amplifiers  can  be  obtained  up  to  stage  gains  of  approximately  100. 

The  effects  of  ASE  are  exacerbated  by  the  inhomogeneously  broadened  C02  laser 
transition  in  which  hole  burning  on  a  few  rotational  lines  leave  spectral  regions  with  high 
gain  on  adjacent  rotational  lines.  ASE  from  these  lines  can  deplete  the  inversion  prior  to 
collisional  mixing  into  the  hole  burned  regions  and  subsequent  stimulated  emission  from 
these  lines. 

The  effects  of  ASE  on  laser  power  extraction  efficiency  and  stage  gain  are  shown  in 
Fig.  5.4.  For  the  EBL  amplifier  configuration,  curve  2  indicates  that  extraction  efficiencies 
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as  high  as  73%  can  be  achieved  at  a  stage  gain  of  1000  in  the  presence  of  ASE  and  a 
1%  specular  reflection  from  the  cavity  windows  or  other  elements  within  the  optical  train. 
This  is  very  close  to  the  extraction  efficiency  of  73.5%  predicted  by  curve  1  which  does 
not  account  for  ASE.  The  relative  insensitivity  of  the  EBL  amplifier  to  ASE  is  a  direct 
result  of  the  beam  expansion  which  (1)  reduces  the  intensity  feedback  from  a  specular 
reflection  along  an  axial  ray  by  a  factor  of  10  and  (2)  provides  a  stage  gain  of  1000  with 
an  integrated  gain  along  a  ray  of  only  100.  The  remaining  factor  of  10  in  the  stage  gain 
is  given  by  the  area  expansion  ratio.  For  a  conventional  single  pass  amplifier  with  the 
configuration  shown  in  Fig.  5.6,  Curve  3  indicates  the  relationship  between  stage  gain  and 
laser  power  extraction  efficiency  for  no  ASE,  and  Curve  4  indicates  this  relationship  in  the 
presence  of  ASE  and  a  1%  specular  reflection  from  the  windows  or  other  elements  of  the 
optical  train.  The  precipitous  decrease  in  laser  power  extraction  efficiency  for  stage  gains 
above  100  in  a  conventional  single  pass  amplifier  in  the  presence  of  ASE  and  a  realistic 
value  of  specular  reflection  is  clearly  shown  in  Curve  4. 

5.6  STEADY-STATE  CONDUCTION  COOLING 


A  critical  requirement  of  a  C02  laser  is  that  of  cooling.  At  a  repetition  rate  of  u  (Hz), 
an  extraction  of  IV  ( J /liter) ,  a  volume  of  V  (liters)  and  an  electrical  efficiency  of  77,  the 
amount  of  heat  Q  which  must  be  carried  away  from  the  laser  to  maintain  it  at  constant 
temperature  is 


vWV 

V 


(5.3) 


For  example,  for  u  =  10  Hz,  W  =  25  J /liter,  V  =  5  liters  and  77  =  0.1,  values  typical  of  a 
1-Atm  C02  laser,  one  obtains  Q  «  10  kW.  Carrying  away  this  much  heat  requires  a  flow 
loop,  which  increases  system  weight  substantially.  In  this  section  the  alternative  approach 
of  conduction  cooling  will  be  addressed  instead  as  an  advanced  concept  to  minimize  the 
weight  of  the  source  laser.  As  shown  in  Fig.  5.S,  the  heat  may  be  conducted  away 
by  the  sidewalls  along  a  direction  perpendicular  to  the  discharge  direction.  Conduction 
cooling  is  a  viable  option  if  the  laser  is  operated  at  a  low  pressure  and  a  correspondingly 
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Figure  5.8:  Geometry  of  laser 
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low  extraction.  For  example,  for  W  =  1  J/liter,  V  =  10  liters,  the  repetition  rate  and 
efficiency  mentioned  above  may  be  maintained  with  conduction  cooling  of  only  1  kW.  The 
sacrifices  of  lower  extraction  and  increased  laser-head  volume  are  reasonable  when  they 
are  linked  to  the  much  more  significant  weight  and  volume  reduction  made  possible  by  the 
abandonment  of  flow  cooling. 

5.6.1  Conduction  Equations 

For  a  repetitively  pulsed  device  the  temperature,  density  and  pressure  all  vary  during 
each  pulse.  A  repetitively  pulsed  CO2  laser  applicable  to  laser  radar  is  expected  to  have 
pulses  of  width  ~  30  ^sec.  These  pulses  are  so  short  that  conduction  within  the  laser 
gas  may  be  ignored  during  each  pulse.  The  equations  coupling  the  temperature,  density 
and  pressure  may  then  be  reduced  to  a  third-order  equation  for  the  density  disturbance 
driven  by  the  pulsed  heating.  The  finite  acoustic  speed  limits  the  response  of  density 
disturbances  of  scale  sizes  greater  than  ~  1  cm  during  each  pulse.  Between  pulses  the 
situation  is  quite  different,  however.  The  pulses  are  far  enough  (~  100  msec)  apart  for 
gaseous  conduction  of  the  heat  to  be  significant.  Since  the  rate  of  heat  conduction  is 
proportional  to  the  temperature  gradient,  small-scale  temperature  gradients  will  tend  to 
conduct  away,  leaving  behind  predominantly  the  largest-seale  gradients.  After  a  large 
number  of  such  pulses,  the  laser  will  reach  a  quasi-steady  state,  in  which  the  temperature 
and  density  are  periodic  in  time  with  the  period  of  the  electrical  pumping.  In  this  limit, 
the  effect  of  each  individual  electrical  pulse  on  the  beam  quality  is  not  expected  to  be 
severe.  However,  the  cumulative  effect  of  a  large  number  of  pulses,  balanced  against  the 
conduction  loss,  is  expected  to  lead  to  steady-  state  density  and  temperature  profiles  which 
significantly  influence  the  wavefront  of  the  laser  beam.  This  steady-state  effect  may  be 
calculated  approximately  by  ignoring  the  pulsed  nature  of  the  excitation  and  treating  it 
as  if  the  heating  is  uniform  in  time,  with  the  same  average  rate  as  that  of  the  actual 
repetitively  pulsed  heating. 

For  heating  which  is  uniform  in  time,  the  equation  for  the  temperature  T,  including 
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the  effect  of  thermal  conduction  of  the  gas,  is 

r)T 

cp— -«V2T=5  (5.4) 

where  cp  is  the  specific  heat  (at  constant  pressure)  per  unit  volume,  k  is  the  thermal 
conductivity  of  the  gas,  and  5  is  the  heating  rate  per  unit  volume.  Eq.  (5.4)  may  be 
written  in  the  alternative  form'9) 

^-DV2T=-  (5.5) 

dt  cp 

where  D  ~  k/cp  is  the  heat  diffusivity,  which  scales  as  the  product  of  the  acoustic  speed  and 
the  mean-free  path  between  collisions.  It  is  thus  inversely  proportional  to  the  gas  density. 
Therefore,  if  the  heating  source  term  S/cp  is  kept  constant  as  the  density  is  decreased  from 
1  atm,  the  effect  of  heat  conduction  increases  significantly.  This  may  also  be  seen  directly 
from  Eq.  (5.4):  if  one  decreases  the  density  and  simultaneously  decreases  5  proportionally, 
then  since  k  is  independent  of  density  the  conduction  term  takes  on  increasing  significance 
as  the  density  is  lowered. 

An  estimate  of  the  efficacy  of  conduction  cooling  for  relevant  laser  conditions  may 
be  made  directly  from  Eq.  (5.4).  Consider  a  laser  channe1  >f  dimensions  4  x  20  cm2, 
with  the  surfaces  separated  by  d.  ~  4  cm  maintained  at  constant  temperature  To-  For 
a  3:2:1  He:NT2:C02  mix  at  400°Iv,  we  have  k  =  1  x  I0~3W/cm°K.  In  this  geometry  the 
cooling  is  dominated  by  the  effect  of  the  walls  spaced  4  cm  apart.  Neglecting  the  other 
walls,  assuming  a  spatially  uniform  heat  deposition,  and  taking  the  steady-state  limit  ( T 
independent  of  time),  one  may  integrate  Eq.  (5.4)  analytically  to  obtain 

AT(x)  =  jfid  -  x)  (5.6) 

2ktj 

where  x  is  the  distance  from  one  of  the  walls.  At  u  =  10  Hz,  IV  =  10~J  J/c:  r/  =  0.1. 

Eq.  (5.4)  predicts  a  steady-state  temperature  rise  of  ~  200  K  halfway  between  the  walls. 

The  1  J/lit  extraction.  10'Z  efficiency  and  10  Hz  repetition  rate  will  be  taken  to  be 
the  baseline  parameters  of  a  CO2  laser  amplifier  to  be  used  for  laser  radar.  A  preliminary 
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analysis  indicates  that  the  d  =  4  cm  spacing  gives  adequate  conduction  cooling  at-  v  =  10 
Hz;  therefore,  unless  otherwise  noted,  this  spacing  will  be  tentatively  considered  a  baseline 
parameter,  as  well.  It  is  anticipated  that  at  higher  repetition  rates  a  smaller  wall  spacing 
will  be  required.  An  alternative,  more  conservative  baseline,  with  d  =  2.4  cm,  will  also  be 
referred  to  in  this  report. 

5.6.2  Thermal  Effect  of  Pulsed  Nature  of  Heating 

One  can  also  solve  Eq.  (5.4)  for  the  temperature  as  a  function  of  time,  assuming 
that  the  energy  deposition  in  the  medium  is  uniform,  and  that  the  changes  take  place  at 
constant  pressure.  One  obtains  a  maximum  temperature  rise  A Tm  given  by*10) 

ATm  =  A Tp  1  +  (-)  (—  _  x  “  3(e9v  _  i)  +  •  •  •) 

In  Eq.  (5.7)  A Tp  =  W jepr]  is  the  temperature  rise  due  to  a  single-pulse  energy  deposition. 


and  the  dimensionless  constant  u  is  given  by 


K  /7T\2 


where  cp  is  the  specific  heat  at  constant  pressure.  For  d  =  4  cm,  one  obtains  v  «  1.5. 
A Tp  *  210  K,  and  A Tm  »  300  K. 

Finally,  it  is  useful  to  determine  the  rate  at  which  temperature  nonuniformities  of 
different  scale  sizes  are  dissipated  by  thermal  conduction  between  pulses.  Let  8To{x)  rep¬ 
resent  the  nonuniform  temperature  rise  at  the  end  of  a  pulse.  Expanding  the  temperature 
in  Eq.  (5.4)  in  spatial  Fourier  components  yields  during  the  time  between  pulses 

dSt(k)  ,  ;  2  c-tv  ,  i  n 

c p  — — —  -f -  nfc  o i  \k)  —  0 


with  the  solution 


6T(k)  =  6t0(k)e-Kk2t/c> 


At  a  repetition  rate  of  v  Hz.  the  pulses  are  separated  by  a  time  t  =  1/.'.  Therefore  the 
residual  temperature  rise  of  spatial  frequency  k„  =  nn/d  due  to  a  single  previous  pulse  is 
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^  =  e--  =  exp(^L) 

6Tnn  cpd2u 

where  v  is  defined  in  Eq.  (5.8),  and  n  =  1,2. . .  is  an  integer  which  ^aracterizes  the  spatial 
frequency  as  a  multiple  of  it/d.  As  discussed  above,  v  ss  1.5  when  d  =  4  cm.  Therefore, 
the  temperature  component  with  the  lowest  spatial  frequency  (n  =  1)  supported  by  the 
channel  decays  by  a  factor  of  ~  exp(1.5)  =  4.5  between  pulses.  The  component  with  the 
next  lowest  spatial  frequency  (n  =  2)  decays  by  a  factor  of  ~  exp(6)  =  400  between  pulses. 
Thus  the  small-scale-  size  nonuniformities  quickly  dissipate  away  between  pulses,  and  to 
an  excellent  approximation  only  the  largest-scale  temperature  nonuniformity  needs  to  be 
taken  into  account. 

5.6.3  Scaling  of  Conduction-Cooled  Laser 

Eq.  (5.6)  may  be  used  to  scale  the  laser  size  as  a  function  of  electrical  power  loading. 
The  steady-state  temperature  Tc  halfway  between  the  walls  may  be  expressed  in  terms  of 
the  electrical  power  loading  P'  —  Wu/rj  as 

^  _  A  m  d2P' 

Tc  —  To  =  AT  =  — —  (5.9) 


Thus  the  temperature  rise  AT  is  proportional  to  the  power  loading  and  the  square  of  the 
distance  between  the  walls.  This  relationship  is  depicted  graphically  in  Fig.  5.9,  which 
shows  the  temperature  rise  as  a  function  of  the  power  loading  for  various  values  of  channel 
width  (wall  separation). 

The  scaling  with  pulse  repetition  frequency  (PRF)  v  may  be  obtained  by  rewriting  Eq. 


(5.9)  in  the  form 


Sac  Tc  —  To 

(ii)~ 


(5.10) 


In  Fig.  5.10  the  relationship  between  channel  width  and  PRF  is  plotted  for  various  values 
of  temperature  rise,  assuming  a  single-pulse  electrical  energy  loading  W/t]  of  10  J /lit.  It 
is  seen  that  at  10  pps  and  a  4-cm  channel  width  the  temperature  rise  is  200°C. 
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Figure  5.10:  Maximum  channel  width  versus  PRF  for 
various  steady-state  temperature  rises 
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I 

I 

Higher  repetition  rates  will  require  smaller  channel  widths  in  order  to  maintain  the 
|  same  temperature  rise.  Table  3  lists  steady-state  temperatures  and  transient  temperature 

rises  for  a  laser  head  with  a  more  conservative  2.4-cm  channel  spacing  operated  at  various 
I  pulse  repetition  rates  (PRF’s).  At  the  highest  PRF  shown  (50  Hz)  it  is  necessary  to 

|  maintain  the  channel  walls  at  a  temperature  of  200°K  (as  indicated  in  Table  3)  in  order 

to  prevent  the  gas  from  reaching  temperatures  at  which  the  lasing  efficiency  would  be 
degraded.  The  200°K  temperature  is  still  well  above  the  boiling  point  (~  160-170°K)  of 
C02  over  the  estimated  pressure  range  (30-80  torr)  which  will  apply  during  operation  of 
I  the  laser. 


Table  3:  Temperature  Variation  for  2.4-cm- Channel  Laser  Head 


Steady  State  Temperature  (K) 

Transient  Temperature  Rise  (K) 

Repetition 

At  Electrode 

At  Channel 

Peak 

Beginning  of 

Rate  (Hz) 

Center 

Pulse 

10 

295 

370 

210 

0 

20 

295 

440 

250 

30 

50 

200 

560 

420 

220 

5.6.4  Impact  of  Conduction  Cooling  on  Beam  and  Spectral  Quality 
Width  in  the  Steady  State 

In  the  preceeding  section  it  was  shown  that  for  a  uniform  heat  deposition,  conduction 
cooling  leads  to  a  steady  state  with  a  quadratic  temperature  variation  across  the  channel 
[see  Eq.  (5.6)].  The  quadratic  variation  of  the  temperature  will,  via  the  equation  of  state. 


lead  to  a  corresponding  variation  in  density  p : 


p  —  A  [  Tq  + 


W  vx{d  —  x) 


2  KT) 


-l 


(5.11) 


where  the  constant  A  may  be  determined  from  the  constraint  that  the  mean  value  of  p, 
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averaged  over  x,  must  be  independent  of  the  heating.  The  main  aberration  coming  from 
this  density  variation  will  be  a  negative  cylindrical  focus  along  x ,  i.e.,  perpendicular  to  the 
channel  walls.  The  associated  radius  of  curvature  R  may  be  calculated  from 


1_ 

R 


(5.12) 


where  A  is  the  wavelength,  <t>  is  the  phase  aberration,  and  V2<£  is  evaluated  at  x  =  d/2, 
halfway  between  the  channel  walls.  The  relationship  between  phase  and  density  is 


2ttLi  /3g 

A  Pa 


(5.13) 


where  f3g  is  the  Gladstone- Dale  constant  (the  refractive  index  minus  unity  at  standard 
temperature  and  pressure  [STP])  of  the  laser  gas,  p ^  is  the  gas  density  at  STP,  and  Lj  is 
the  propagation  length.  Combining  Eqs.  (5.11)-(5.13)  yields 


2 ~  —  3  p0  ^ 

R  ~  3  d*P3pAl+p. 


(5.14) 


where  p0  is  the  initial  density  of  the  gas,  and  p  is  given  by 


Tc  -  r0  _  Wud} 
Tq  8k7]Tq 


(5.15) 


Here  To  is  the  wall  temperature  and  Tc  is  the  temperature  at  the  center.  The  radius  of 
curvature  for  a  single  pass  through  the  amplifier  ( L\  =  125  cm)  when  po/ Pa  =  0.04  and 
p  =  1  is 

R  %  120  m 


[  which  is  reasonably  large.  For  4  passes  through  the  amplifier,  with  no  correction  between 

; 

passes,  the  radius  of  curvature  would  still  be  ~  30  m,  corresponding  to  ~  0.03  Diopters. 
The  cylindrical  focus  error  can  easily  be  compensated  statically  by  giving  the  laser  beam 
incident  on  the  laser  amplifier  a  cylindrical  focus  along  x.  If  the  amplifier  utilizes  an 
expanding  beam  geometry  then  the  optics  producing  the  expansion  (perpendicular  to  x) 
will  be  in  place,  and  all  that  will  be  needed  is  a  small  modification  of  these  optics  to  produce 
!  5-23 
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the  relatively  small  focus  along  x.  Thus,  the  cylindrical  focus,  which  is  the  dominant 
aberration  due  to  the  steady-state  thermal  loading,  can  easily  be  compensated  and  should 
produce  no  appreciable  beam-quality  degradation. 

Higher-order  aberrations  due  to  the  thermal  loading  may  also  be  present  if  the  electrical 
heating  is  not  uniform  along  x.  These  may  be  estimated  by  taking  the  spatial  Fourier 
transform  of  Eq.  (5.4)  in  the  steady  state,  yielding 

T(k)  =  SS  (5.16) 

The  amplitude  of  the  component  of  the  temperature  of  spatial  frequency  k  is  thus  shown  to 
be  proportional  to  the  heating  nonuniformity  of  the  same  spatial  frequency  and  inversely 
proportional  to  the  square  of  the  spatial  frequency.  Therefore,  as  discussed  above,  small- 
scale  nonuniformities  of  the  steady-state  heating  do  not  produce  strong  temperature  and 
density  nonuniformities.  The  reason  is  that  for  a  given  temperature  difference,  heat  is 
conducted  more  rapidly  the  smaller  the  scale  size  of  the  temperature  difference.  In  the 
steady  state  the  dominant  density  disturbance  will  always  be  the  low-order  cylindrical 
focus,  which,  as  discussed  above,  is  easily  compensated.  If  it  turns  out  to  be  necessary, 
it  is  also  possible  to  compensate  some  of  the  higher-order  aberrations  by  utilizing  a  static 
adaptive  optic  which  may  be  adjusted  in  place  to  minimize  the  laser  divergence. 

It  should  be  noted  that  the  steady-state  temperature  and  density  profile  distributions 
have  no  effect  on  the  spectral  purity  of  the  laser.  Only  the  pulsation  of  the  heating,  which 
causes  time-dependent  changes  in  the  phase,  can  broaden  the  spectrum  by  producing  chirp 
and  frequency  noise. 

A  preliminary  analysis,  sketched  above,  thus  indicates  that,  as  far  as  the  steady  state 
is  concerned,  conduction  cooling  to  channel  walls  is  a  viable  technique  for  carrying  away 
the  heat  generated  in  a  pulsed  CO2  laser  amplifier  operated  at  a  low  pressure  (0.04  atm). 
The  heat  generated  by  each  electrical  pulse  is  sufficient  to  raise  the  temperature  of  the 
laser  gas  significantly,  on  the  order  of  ~  300  °C.  However,  most  of  this  heat  conducts  away 
to  the  walls  between  pulses.  After  several  laser  pulses,  a  quasi-steady-state  temperature 
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and  density  profile  appears.  It  is  comparable  with  the  profile  that  one  would  obtain  if  the 
pulses  were  replaced  by  continuous  heating  at  the  same  rate.  The  profile  is  dominated  by 
a  quadratic  cylindrical  focus  which,  for  four  passes  through  the  amplifier,  is  expected  to 
be  ~  0.03  Diopters.  This  negative  focus  is  transverse  to  the  beam-expansion  direction, 
and  is  easily  removed  with  a  slight  modification  of  the  focussing  optics.  It  is  noteworthy 
that  the  negative  focus  does  not  appear  until  the  laser  amplifier  is  operated  for  several 
pulses  so  that  it  reaches  its  steady-state  thermal  balance.  Higher-order  aberrations  due  to 
the  steady-state  heating  are  predicted  to  be  small,  since  the  temperature  gradients  which 
drive  them  are  strongly  dissipated  by  heat  conduction.  However,  it  is  possible  to  remove 
these  aberrations,  as  well,  by  utilizing  a  relatively  simple  static  adaptive-optics  system. 
The  steady-state  heating  should  have  no  effect  on  the  spectral  purity  of  the  amplified  laser 
beam. 

5.7  ELECTRICALLY  GENERATED  DENSITY  DISTURBANCES 

The  electrical  energy  used  to  pump  the  laser  medium  will  unavoidably  have  some  degree 
of  spatial  nonuniformity.  This  will  result  in  the  generation  of  acoustical  and  entropy 
disturbances,  which  will  in  turn  produce  refractive-index  nonuniformities,  an  increased 
spectral  width,  and  beam-quality  degradation  in  the  laser.  Because  the  channel  walls 
of  the  discharge-pumped  laser  are  solid,  there  is  no  intrinsic  spatial  discontinuity  in  the 
pump-energy  distribution,  such  as  that  which  occurs  in  an  e-beam-pumped  laser.*11)  Still, 
if  the  laser  medium  of  interest,  discharge-pumped  by  ~  10  J /lit  in  a  ~  30-psec  pulse,  were 
used  as  an  oscillator,  it  would  be  extremely  difficult  to  maintain  the  pumping  uniform 
enough  to  make  the  oscillator  beam  quality  and  spectral  quality  acceptable.  However,  as 
will  be  shown,  the  beam-quality  degradation  and  frequency  noise  are  not  so  severe  if  the 
laser  medium  is  used  instead  as  an  amplifier.  It  is  expected  that  the  dominant  source 
of  the  discharge  nonuniformities  will  be  nonuniform  UV  preionization.  It  is  thus  possible 
to  make  the  laser  light  average  over  the  discharge  nonuniformity  by  ensuring  that  the 
UV-preionizing  structure  changes  along  the  direction  of  propagation  of  the  laser  light. 
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The  distance  over  which  this  structure  must  change  significantly  in  order  to  ensure  good 
beam  quality  and  frequency  purity  will  determined  in  detail.  In  this  report,  a  preliminary 
estimate  of  the  requirements  for  good  beam  quality  will  be  made. 

5.7.1  Density  Response  to  Nonuniform  Pumping 

The  geometry  of  the  laser  is  shown  in  Fig.  5.S.  The  dominant  direction  of  discharge 
nonuniformity  is  expected  to  lie  along  the  direction  x  transverse  to  the  channel  surfaces. 
Thus,  the  contributions  of  nonuniform  pumping  to  density  variations  along  the  discharge 
direction  y  and  the  optic  axis  z  may  be  neglected,  and  to  a  good  approximation  the 
perturbation  p  of  the  density  satisfies  a  one-dimensional  wave  equation  of  the  form(12* 


d_ 

dt 


(5.17) 


where  Q(x,t)  is  the  heating  rate,  u  is  the  speed  of  sound,  and  7  is  the  ratio  of  specific 
heats  of  the  lasing  gas.  To  determine  the  relevant  characteristics  of  the  solutions  to  these 
equations  one  may  impose  boundary  conditions  p  =  0,  Q  —  0  at  the  walls  x  =  0,  x  —  d, 
which  imply  that  p  and  Q  may  be  expanded  in  discrete  Fourier  series  as 


p{x,t)  =  ^2  pn(t)  sm(knx) 

n 

Q(x,t)  =  Y2Qn(t)sm(knx) 


(5. IS) 


where  k„  =  rcn/d  (n  =  1,2,...)  are  the  allowed  spatial  frequencies.  These  are  not  the 
most  general  boundary  conditions,  but  the  pertinent  characteristics  of  the  solutions  will 
not  depend  on  the  boundary  conditions  except  as  discussed  specifically  below. 
Substituting  Eq.  (5.18)  into  Eq.  (5.17)  yields 


d_ 

dt 


+  (M2  Pn  =  -(7  “  l)KQn{t) 


(5.19) 


A  reasonable  approximation  is  to  consider  the  heating  rate  to  be  constant  over  the  interval 
0  <  t  <  t,  where  r  is  the  pulsewidth  of  the  electrical  pumping.  Thus,  Q(x,t)  may  be 
expressed  as  H(x)/r ,  where  H(x)  is  the  total  electrical  single-pulse  energy  loading  per 
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unit  volume.  Similarly,  Qn(t)  =  Hn/T.  The  solution  of  Eq.  (5.19)  during  the  pulse  is  then 

~uknt  —  sin  (uknt)' 


p,(t)  =  -  (  ]  Hk, 


U‘ 


uknT 


(5.20) 


where  H  is  the  mean  energy  loading  per  unit  volume,  given  by 


(5.21) 


and 

hn  =  ^  J  sin(k„x)dx  (5.22) 

is  the  normalized  nth  Fourier  component  of  the  energy  loading  distribution. 

It  is  useful  to  consider  the  form  of  Eq.  (5.20)  in  the  short-time  and  long-time  regimes. 
kn  =  7r njd  is  the  frequency  of  the  ntk  acoustic  mode  of  the  channel.  For  each  mode  n  of 
wavelength  2ir/kn  =  2 d/n,  one  may  make  the  following  two  observations:  (i)  When  not 
enough  time  has  passed  for  a  sound  wave  to  have  travelled  across  an  acoustic  wavelength 
(uknt  <  1),  sin uknt  %  uknt  so  that  p„  ~  t3.  This  is  the  well-known  transient  “f3”  regime 
of  thermal  blooming,  (ii)  For  longer  times  (uknt  >>  1),  corresponding  to  many  acoustic 
transits  across  an  acoustic  wavelength,  the  blooming  reaches  the  steady  state,  and  pn  is 
instead  proportional  to  t.  Thus,  the  density  will  respond  more  strongly  to  small-scale 
energy-loading  nonuniformities  (with  ufc„r  >>  1)  than  to  larger-scale  nonunifc.mities 
(with  uknr  <  1). 

In  the  following,  spatial  averages  (over  z)  will  be  denoted  by  <  ...  >x,  and  time 
averages  by  <  . . .  >t.  As  will  be  shown  in  the  next  section,  the  beam  quality  is  determined 
by  the  spatial  variance  of  the  density  perturbation,  averaged  over  the  pulse.  Using  Eq. 
(5.20),  the  time-averaged  variance  may  be  expanded  into  its  spatial  Fourier  components, 
with  the  result 


«p 2  >*><= 


where  B ,  which  has  the  dimensions  of  a  density,  is  given  by 

'(7-1)  H' 


B  = 


(5.23) 


(5.24) 
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and  /„  =  1(9  =  uknr )  is  the  response  function  of  the  nth  acoustic  mode  of  the  channel  to 
the  energy-  loading  pulse  of  length  r: 


Q 3  1  1 

— —  2  sin  9  4-  -9(1  +  9  cos  9)  —  —  sin  29 


-  (5-25) 


The  behavior  of  I  as  a  function  of  the  normalized  pulsewidth  9  is  shown  in  Fig.  5.11.  It  is 
seen  that  the  transition  between  the  short-  and  long-time  regimes  is  relatively  sharp,  and 
that  a  good  approximation  for  the  function  1(9)  is 


f  9s-/ 252,  if  9  <  3; 
\b  * >  3. 


(5.26) 


The  f4  dependence  at  short  times  corresponds  to  the  t6  dependence  of  p2  implicit  in  Eq. 
(5.20),  corrected  for  the  fact  that  at  a  constant  heating  rate  B 2  would  be  proportional  to 
t2. 


5.7.2  Beam  Quality 

The  energy  loading  nonuniformity  and  the  density  disturbance  which  it  causes  will  in 
general  not  be  independent  of  position  z  along  the  optic  axis  of  the  laser.  In  fact,  the 
nonuniformities  will  occur  in  layers  along  z,  with  successive  layers  spatially  uncorrelated 
with  each  other.  The  phase  aberrations  from  each  layer  will  therefore  add  randomly  to 
one  another,  much  like  a  random  walk.  The  following  calculation  will  concentrate  first  on 
the  effect  of  a  single  layer.  The  addition  of  many  layers  will  be  carried  out  subsequently. 

The  phase  aberration  <p(x,t)  seen  by  a  ray  propagating  along  z  across  a  single  layer  is 
directly  proportional  to  p(x,<): 


d(-M)  =  ap(x,<)  =  z^-I^-p(x,t)  (5.27) 

A  Pa 

where  Lq  is  the  correlation  length  along  z  (i.e.,  the  thickness  of  the  layer),  A  is  the  wave¬ 
length  of  the  laser  light,  /3g  is  the  Gladstone-Dale  constant  (the  difference  between  the 
refractive  index  of  the  amplifying  medium  and  unity  at  standard  temperature  and  pres¬ 
sure  [STP]  ),  and  pA  is  the  density  of  the  amplifying  medium  at  STP.  6(x.  t)  may  be 
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expanded  in  a  Fourier  series  of  the  form 


<p(x,t)  =  Y^^n(t)sinknx 


(5.2S) 


where  (pn  —  apn ,  and  the  spatial  variance  of  the  phase  aberration,  averaged  over  time,  may 
similarly  be  expanded  as 

<<  4>2  >z>t=  X]  hnIn  (5'2i)) 

^  n 

Here  $  =  cxB  is  given  explicitly  by  the  following  convenient  expression: 


,  2  ttL0 


?-n  (h_ 

7  )  \pa 


(5.30) 


where  pA  is  the  standard  pressure  (1.01  x  105  Nt/m2,  or,  in  energy  density  units,  101 
J /lit.)  Note  that  $  depends  only  on  the  energy  loading  per  unit  volume  and  the  cell 
length;  in  particular,  it  is  independent  of  pressure.  However,  at  low  gas  pressures  and 
correspondingly  low  energy  loading,  <£  can  be  reasonably  low.  A  typical  value  of  $  is  0.31 
radians  at  £0  =  5  cm,  A  =  10  pm,  /3g  =  3  x  10-4,  7  =  1.5,  and  H  =  10  J /lit. 

Beam  quality  is  determined  by  the  relative  value  of  the  peak  far-field  intensity,  known 
as  the  Strehl  ratio. In  principle,  the  Strehl  ratio  depends  on  both  the  phase  and  intensity 
nonuniformities  in  the  near  field.  In  practice,  however,  the  contribution  of  the  intensity 
nonuniformity  to  deviations  of  the  Strehl  ratio  from  unity  is  usually  negligible.  For  the 
case  at  hand  the  Strehl  ratio  S(t)  may  then  be  written  as 

1 2 


S(t)  = 


<  eWxA)  >r 


(5.31) 


where,  as  above,  the  brackets  <  ...  >x  denote  a  spatial  average  over  the  near-field  beam 
(i.e.,  over  x).  Expanding  the  exponent  in  powers  of  <t>  yields 

5(t)  «  1  —  (<  <t>2  >z  -  <  <t>  >1)  (5.32) 

to  order  <f>2.  In  general,  one  is  interested  in  the  temporal  average  of  the  beam  quality.  Let 
5  be  the  time  average  of  S(t).  It  follows  from  Eq.  (5.32)  that 

5  %  1-  <<  (j>2  >z>t  +  «  <$>  >l>i  (5.33) 
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When  \<j>—  <  <f>  >  |  >  1  this  approximation  breaks  down,  and  higher  powers  of  -<j>  must 
be  included.  If  the  energy-loading  distribution  satisfies  Gaussian-normal  statistics  then 


the  power  series  may  be  summed,  yielding 


S  «  exp  [-  <<  <j> 2  >i >t  +  <<  <t>  >*><] 


(5.34) 


[For  other  statistical  distributions  of  the  energy  loading  Eq.  (5.34)  is  correct  only 
to  order  <p 2;  Eq.  (5.34)  and  all  subsequent  equations  describing  beam  quality  should  in 
general  be  considered  to  be  correct  only  to  this  order.] 

The  presence  of  hard  boundaries  confining  the  gas  in  a  fixed  volume  forces  the  mean 
density,  averaged  over  x,  to  be  a  constant,  independent  of  t.  Therefore  <  <f>  >x=  0. 
Combining  Eqs.  (5.26),  (5.29),  and  (5.30),  one  arrives  at  the  following  result: 


i.  .  R<no  n>no 


where 


n°  7T  UT 


(5.35) 


(5.36) 


is  approximately  the  number  of  acoustic  transits  across  the  channel  during  the  electrical 
pulse.  By  way  of  example,  for  d  =  4  cm,  r  =  30  ^sec,  and  u  ss  400  m/sec,  the  value  of  no 
is  ~  3.2. 

It  should  be  noted  that  beam  quality  is  often  measured  in  “number  of  times  diffraction 


limited”  ( XDL ).  The  relationship  between  XDL  and  5  is 


XDL  = 


(5.37) 


so  that  Eq.  (5.35)  implies  that 


(  L  n<no  n>no 


(5.3S) 


The  above  results  apply  to  the  propagation  through  a  single  correlation  length  (along 
r)  of  the  density  disturbance.  Since  successive  layers  are  spatially  uncorrelated,  the  prop¬ 
agation  through  several  correlation  lengths  results  in  the  addition  of  the  squares  of  the 
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phase  errors.  The  generalizations  of  Eqs.  (5.35)  and  (5.38)  are  thus 

5  *  «p  e  n'h» + 5  e  *’.] }  (5M> 

k  L  n<i»o  n>no  J  ) 

and 

A'£>L*exp  j  +  {^pj  +  (5.40) 

l  ®  .  n<no  n>no  .  ) 

where  L  is  the  total  distance  travelled  by  the  laser  beam  through  the  amplifier,  so  that 
L/Lq  is  the  number  of  correlation  lengths  traversed  by  the  laser  beam. 

It  is  useful  to  break  the  right-hand  side  of  Eq.  (5.40)  into  two  parts,  as  follows: 


(5.40) 


XDL  as  exp 


[8l  +  <5s] 


(5.41) 


where  81  and  6s  are  the  contributions  from  large-  and  small-scale  nonuniformities,  respec¬ 


tively,  in  the  electrical  loading: 


*  -  (=r)  £ 

n<n0 


(5.42) 


81  is  the  contribution  of  transient  thermal  blooming  caused  by  large-scale  nonuniformities 
(scale  size  >  ut)  in  the  electrical  energy  loading.  Because  of  the  factor  of  ( ~uTn/d)i 
in  this  term,  it  is  generally  small  (see  Fig.  5.10).  On  the  other  hand,  8s  is  due  to  the 
steady-state  thermal  blooming  caused  by  small-scale  nonuniformities  (scale  size  <  ut)  in 
the  energy  loading.  Because  of  the  relatively  strong  steady-state  response  of  the  medium 
to  the  smaller  scales,  it  is  65  which  must  be  carefully  monitored  in  the  design  of  the  laser 
amplifier. 


A  useful  alternative  expression  for  85  is 

*S  *  |  <  hi  >x=  |[/l5,rm,]2 

where  hs(x)  is  the  small-scale  fractional  modulation  in  the  energy  loading: 


(5.43) 


hs(x)  = 


H(x)-H 


(5.44) 
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Combining  Eq.  (5.44)  with  the  definition  of  $  yields 

.Y£L«exp  Ho[/is,rm,j2|  (5.45) 

For  a  given  required  XDL,  Eq.  (5.45)  may  be  used  to  determine  the  acceptable  range  of 
Lq  and  rms  variation  in  the  energy  loading.  By  way  of  example,  for  three  passes  through 
a  125-cm  amplifier  one  has  L  =  375  cm.  With  A  =  10  fim,  0g  =  2  x  10-4,  7  =  1.5,  and  for 
H  =  10  J/liters,  Eq.  (5.45)  reduces  to 

XDL  »  exp  |  +  Xl(-^9rmj)2-|  (5.46) 

where  Lq  is  in  cm.  Eq.  (5.46)  implies  that  energy- loading  nonuniformities  as- large  as 
~  20  —  30%  rms  are  consistent  with  excellent  beam  quality  (1.1  XDL)  if  the  correlation 
length  is  of  the  order  of  10  cm  or  less. 

This  preliminary  analysis  thus  indicates  that  the  requirements  for  little  effect  on  beam 
quality  from  the  acoustics  driven  by  a  single  pulse  are  relatively  easy  to  achieve.  A  more 
detailed  analysis,  including  realistic  deviations  from  Gaussian- normal  statistics,  will  be 
carried  out  during. 

5.7.3  Frequency  Noise 

A  second  important  effect  of  the  density  changes  during  the  pulse  is  changes  in  the 
frequency  of  the  amplified  laser  light.  In  general,  these  frequency  changes  may  be  divided 
into  two  distinct  parts:  chirp  and  frequency  “noise.”  The  chirp  i/c*  is  the  change  in  the 
instantaneous  central  frequency  of  the  light,  while  the  noise  Su  is  associated  with  deviations 
from  the  instantaneous  central  frequency,  i.e.,  creation  of  Fourier  components  at  other 
frequencies.  For  a  laser  amplifier  with  solid  walls,  the  chirp  should  be  nearly  zero  if  the 
laser  beam  uniformly  fills  the  entire  region  between  the  walls  (i.e.,  if  the  laser  intensity  is 
independent  of  x).  Deviations  from  uniformity  may  produce  some  chirp. 

5.7.4  Mode-Medium  Instability 

Another  effect  that  could  in  principle  prove  relevant  to  controlling  the  beam  quality  of  a 
pulsed  CO;  laser  is  the  mode-medium  instability  (MMI),  which  has  been  described  in  detail 
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in  Ref.  11.  The  MMI  is  due  to  the  interaction  between  the  laser  intensity  distribution,  the 
dependence  of  the  heating  rate  on  the  laser  inversion,  the  acousto-density  response  of  the 
gaseous  medium  to  the  heating  distribution,  and  the  refraction/diffraction  of  light  by  the 
density  disturbance. 

It  is  anticipated  that  a  CO2  laser  amplifier  of  the  type  used  for  laser  radar  will  not 
suffer  any  significant  loss  from  MMI.  As  was  shown  by  Flusberg  et.  al.,  the  instability 
is  important  in  a  laser  oscillator/12'  but  does  not  cause  a  loss  of  beam  quality  in  a  laser 
amplifier/14,15'  In  addition,  even  in  an  oscillator,  the  rate  at  which  the  instability  grows 
depends  on  density;  at  a  density  of  30-40  torr  the  degradation  in  oscillator  beam  quality 
should  occur  on  a  time  scale  of  ~  100  ^isec,  hence  there  should  be  no  significant  effect  on 
the  beam  quality  of  a  30-nsec  pulse. 

5.7.5  Persistent  Acoustic  Disturbances 

Another  important  effect  on  beam  quality  is  that  of  persistent  acoustic  disturbances, 
which  are  created  by  the  heat  deposition  during  each  pulse.  Between  pulses  these  distur¬ 
bances  relax  from  the  combined  effects  of  viscous  damping  in  the  gas  and  imperfect  acoustic 
reflections  from  the  solid  boundaries.  For  a  repetitively  pulsed  laser,  each  successive  pulse 
“feeds”  the  various  acoustic  modes  of  the  laser  cavity.  Each  mode  decays  between  pulses 
but  is  driven  again  by  successive  pulses.  The  result  is  a  steady  state  between  the  driving 
impulses  due  to  the  pulsed  heat  deposition  on  the  one  hand  and  the  damping  by  viscosity 
and  finite  reflectivity  on  the  other.  For  a  small  enough  wall  reflectivity  the  acoustics  should 
damp  away  so  rapidly  that  no  significant  buildup  should  take  place. 

The  effect  of  persistent  acoustic  disturbances  will  be  analyzed,  and  the  wall  reflectivity 
loss  required  to  limit  the  effect  on  beam  quality  will  be  determined. 

5.8  ALL  SOLID  STATE  RF  C02  LASER  PULSER 

RF  pulse  generators  with  output  powers  of  several  tens  of  megawatts  are  not  presently 
available  in  compact  lightweight  packages  (i.e.  <  30  kg).  SRL  proposes  employing  magnetic 
switches  which  will  result  in  a  compact,  lightweight  RF  pulse  generator  capable  of  driving 
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the  laser  head  discussed  in  the  preceeding  sections  of  this  chapter. 

The  use  of  COj  lasers  in  laser  radar  has  sparked  new  interest  in  improving  the  reliabil¬ 
ity  of  discharge  lasers  and  their  associated  pulsed  power  drivers.  Early  attempts  at  driving 
the  lasers  directly  with  thyratrons  or  spark  gaps  have  resulted  in  modulator  lifetimes  of 
only  a  few  hundred  hours.  All-solid-  state  nonlinear  magnetic  drivers  have  greatly  allevi¬ 
ated  this  problem.  The  use  of  nonlinear  magnetic  pulse  compressor  technology  has  allowed 
operation  at  repetition  rates  up  to  several  kilohertz  with  greatly  improved  reliability  and 
reduced  maintenance.  SCR-commutators  combined  with  nonlinear  magnetic  pulser  com¬ 
pressors  allow  for  the  development  of  an  all  solid  state  pulse  generator.  The  basic  principle 
underlying  nonlinear  magnetic  driver  operation  involves  using  a  saturable  core  which  is  an 
inductor  in  a  resonant  circuit.  In  the  unsaturated  state  the  saturable  inductor  exhibits 
a  large  inductance  and  once  saturated  the  inductance  is  reduced  by  the  change  in  per¬ 
meability  of  the  saturable  core.  This  change  in  inductance  is  typically  a  factor  of  several 
thousand.  The  core  is  designed  to  saturate  at  the  peak  of  the  charging  current  waveform 
from  the  previous  stage  and  once  saturated  rapidly  transfers  the  stored  energy  to  the  fol¬ 
lowing  stage.  In  so  doing  the  energy  is  conserved  but  temporally  compressed  resulting  in 
increased  peak  power.  This  nonlinear  saturation  phenomenon  shifts  the  frequency  of  this 
resonant  circuit  by  the  square  root  of  the  permeability  shift  as  the  core  saturates.  These 
stages  are  typically  cascaded  (Fig.  5.12)  and  energy  is  coupled  faster  and  faster  from  one 
stage  to  the  next. 

A  possible  RF  pulser  design  would  be  a  frozen  wave  line  to  build  an  all  solid  state  pulse 
power.  A  frozen  wave  line  consists  of  a  transmission  line  distributed  into  many  isolated 
segments.  These  segments  are  individually  charged  to  potentials  corresponding  to  those 
of  the  desired  output  wave.  The  segments  are  then  simultaneously  connected  together 
allowing  the  wave  to  propagate  along  the  line  to  the  laser  load. 

5.8.1  Conceptual  Design  of  an  All  Solid  State  Pulser 

The  use  of  CO2  lasers  in  laser  radar  has  sparked  new  interest  in  improving  the  reliability 
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Figure  5.12:  Typical  magnetic  switch  operation 
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of  CO2  lasers  and  their  associated  pulsed  power  drivers.  Early  attempts  at  driving  the 
lasers  directly  with  thyratrons  resulted  in  significant  energy  dissipation  in  the  thyratron 
which  led  to  thyratron  lifetimes  of  only  a  few  hundred  hours. 

Nonlinear  magnetic  pulse  compressors  have  greatly  alleviated  this  problem.  The  output 
stages  of  these  units  are  sufficiently  rugged  to  dissipate  the  required  fraction  of  the  drive 
energy  reflected  from  the  nonlinear  load.  However,  thyratrons  are  still  used  in  many  cases 
to  supply  the  input  commutation  for  these  nonlinear  magnetic  pulse  compressors  and 
represent  a  weak  link.  The  use  of  nonlinear  magnetic  pulse  compressor  technology  has 
greatly  eased  the  dl/dt  problems  for  the  thyratrons  but  the  high  average  current  still  leads 
to  a  thyratron  tube  life  which  is  unacceptable  with  respect  to  reliability  and  maintenance. 

The  use  of  solid  state  devices  in  place  of  thyratrons  and  spark  gaps  is  the  key  to 
developing  a  reliable  pulse  power  system.  The  nonlinear  magnetic  circuits  which  pro¬ 
vide  the  necessary  temporal  compression  are  inherently  reliable  as  they  contain  no  active 
elements.  Both  high  pressure  (spark  gaps)  and  low  pressure  (thyratrons)  gas  discharge 
devices  have  a  finite  lifetime.  Electrode  erosion  with  subsequent  gas  contamination  and 
insulator  degradation  will  eventually  degrade  the  performance  of  these  devices.  SCRs  like 
nonlinear  magnetic  compressors  do  not  rely  on  plasma  formation  or  material  displacement 
to  function.  An  SCR  operating  within  its  ratings  will  essentially  last  forever.  Solid  state 
devices  were  originally  developed  by  Bell  Laboratories  for  space  based  applications  where 
maintenance  was  not  possible.  Solid  state  switches,  however,  will  only  operate  outside 
their  ratings  once.  They  are  not  very  forgiving  in  this  respect,  for  exceeding  the  maxi¬ 
mum  voltage,  rms  current  or  dl/dt  ratings  even  instantaneously  will  result  in  irreversible 
damage. 

In  this  effort,  Science  Research  Laboratory  plans  to  exploit  recent  progress  in  Silicon 
Controlled  Rectifier  (SCR)  and  saturable  core  magnetic  switch  technology  to  develop  all 
solid  state  Rr  drivers  which  meet  all  of  the  performance  specifications  for  advanced,  com¬ 
pact  lightweight  CO2  laser  radar  systems.  SCRs  can  be  used  in  conjunction  with  metglass 
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saturable  core  magnetic  switches  to  develop  all  solid  state  modules.  In  addition  to  the 
substantial  SCR  size  advantage  over  thyratrons,  SCRs  do  not  require  filament  transform¬ 
ers,  bias  supplies  and  considerable  other  ancillary  equipment.  The  SCRs  size  advantage 
is  predominately  derived  from  the  increased  density  of  conduction  electrons  in  silicon  as 
opposed  to  the  low  pressure  gas  in  the  thyratron.  In  addition  these  all-solid- state  SCR- 
switched  drivers  can  be  engineered  to  operate  for  greater  than  10n  shots  at  greater  than 
95%  efficiency  to  minimize  waste  heat  rejection. 

The  design  of  efficient  high  gain  pulse  compression  stages  relies  on  the  careful  control 
of  the  magnetic  fields  surrounding  the  core  volume.  Optimization  of  this  packing  factor 
is  crucial  in  magnetic  switch  design  and  is  accomplished  by  enclosing  the  core  in  a  tightly 
fitting  conducting  housing.  The  multiple  turns  axe  formed  as  coaxial  transmission  lines 
which  pass  through  this  housing.  Circulating  currents  set  up  in  this  housing  exclude  the 
magnetic  flux  and  contain  it  in  the  desired  volumes. 

The  design  of  efficient  high  gain  pulse  compression  stages  relies  on  the  careful  control 
of  the  magnetic  fields  surrounding  the  core  volume.  The  losses  in  a  compression  stage  are 
linearly  proportional  to  the  volume  of  core  material  used  and  it  can  be  shown  that  the 


core  volume  requirement  for  this  saturable  inductor  is  given  by: 


volume  =  gain2  •  Pulse  Energy 


Mo  •  7T 


(5.47) 


~  \4  •  (ABS  •  pf)2/  v  ' 

where  A B3  is  the  available  flux  swing  of  the  material  (with  appropriate  biasing  given  by 
+BS  =  2 B3),  gain  is  defined  as  the  temporal  compression  factor  (rciiargc/Tdi,(:hargi), 

and  pf  is  the  packing  factor.  The  packing  factor  is  defined  as 

Sy, 

j„—  hhZ  (5'4S) 

Jail  space 

where  V f  is  the  actual  volume  occupied  by  the  ferri-  (ferro)  magnetic  core  material  exclud¬ 
ing  all  interlaminar  insulation  and  voids.  Optimization  of  this  packing  factor  is  crucial  in 
magnetic  switch  design  and  is  accomplished  by  enclosing  the  core  in  a  tightly  fitting  con¬ 
ducting  housing.  The  multiple  turns  are  formed  as  coaxial  transmission  lines  which  pass 
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through  this  housing.  Circulating  currents  set  up  in  this  housing  exclude  the  magnetic 
flux  and  contain  it  in  the  desired  volumes. 

With  the  saturable  material  completely  enclosed  in  conducting  coaxial  housings,  the 
voltage  per  turn  associated  with  the  dB/dt  of  the  core  appears  across  a  single  narrow  gap 
on  the  outside  diameter  of  each  housing.  Connecting  the  coaxial  transmission  lines  in  the 
inner  and  outer  housings  with  radial  transmission  lines  across  the  top  and  bottom  of  the 
core  completes  the  turns.  The  impedance  of  these  transmission  lines  are  adjusted  so  that 


~  /2  •  La at 

Z LINE  ^  y  g 


(5.49) 


while  simultaneously  every  effort  is  made  to  minimize  Lsat  •  C  which  is  proportional  to 
the  square  of  the  pulse  length.  Here  Lsat  and  C  refer  to  the  saturated  inductance  of  the 
reactor  and  the  value  of  the  storage  capacitance  respectively. 

In  addition  to  optimizing  the  core  structure,  increased  SCR  protection  and  efficiency 
is  provided  by  grading  the  storage  capacitors  throughout  the  chain  so  that  they  gradually 
increase  in  capacitance  in  the  direction  of  the  output.  This  introduces  a  controlled  level  of 
voltage  reversal  per  stage  which  acts  to  slow  the  progress  of  the  reflected  pulse  by  providing 
some  reset.  A  reversed  waveform  cannot  propagate  efficiently  back  through  a  fully-reset 
pulser  as  long  as  the  gain  per  stage  exceeds  a  factor  of  two.  Grading  the  capacitors  in  this 
manner  also  increases  the  efficiency.  Core  losses  introduce  a  resistive  component  to  the 
impedance  which  couples  storage  capacitors  in  each  stage.  Unless  each  successive  capacitor 
is  larger  in  value  than  its  predecessor,  this  resistive  term  will  cause  positive  charge  to  be 
left  on  the  preceding  capacitor  after  the  discharge  cycle  is  completed  and  will  enhance 
propagation  of  pulses  in  the  reverse  direction. 

SRL  is  herein  proposing  to  apply  this  same  technology  to  the  challenge  of  generating 
RF  pulses  in  a  compact  lightweight  driver.  This  RF  generator  will  be  capable  of  supplying 
a  75  joule,  20  MHz  RF  pulse  with  a  30  psec  pulse  duration  at  a  repetition  rate  exceeding 
100  pps. 
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The  basic  concept  of  the  frozen  wave  line  was  originally  conceived  over  50  years  ago 
and  is  illustrated  in  Fig.  5.13.  A  pulse  forming  network  is  divided  into  multiple  segments 
isolated  by  switches.  The  elements  are  charged  from  a  DC  supply  in  such  a  way  that 
the  desired  output  wave  appears  on  the  line.  This  wave  is  in  essence  frozen  into  the  line 
and  prevented  from  propagating  down  the  line  by  the  presence  of  the  open  switches.  The 
switches  are  then  closed  simultaneously  and  the  wave  is  propagates  down  the  transmisssion 
line  to  the  output  terminals.  While  this  idea  appears  simple  enough,  the  switch  technology 
has  not  existed  until  recently  which  could  extend  this  concept  to  high  power  (megawatts) 
at  high  repetition  rate  (100  pps)  in  the  frequency  range  of  interest  (10-20  MHz). 

Recent  developments  in  the  field  of  nonlinear  magnetic  pulse  compression  technology 
provide  a  solution  to  the  switching  problem.  Nonlinear  magnetic  pulse  compression  tech¬ 
nology  is  now  employed  to  drive  CO2  and  rare  gas  halide  excimer  lasers,  copper  vapor 
lasers,  and  free  electron  lasers.  These  modulators  have  all  generated  DC  pulses,  but  now 
we  have  identified  two  techniques  whereby  this  technology  can  be  employed  to  generate 
RF  pulses. 

The  concept  shown  schematically  in  Fig.  5.14  uses  two  separate  nonlinear,  magne¬ 
tically-switched  charging  networks  which  charge  the  isolated  line  elements  to  alternate 
polarities.  Here  the  PFN  network  elements  can  be  segments  of  transmission  lines  but  the 
concept  still  applies  in  the  case  where  simple  capacitors  replace  these  elements  as  shown. 
Transmission  line  segments  have  some  advantages  if  a  square  wave  output  is  desired  as  the 
individual  discrete  line  segments  can  be  center  charged  rather  than  end  charged.  Saturable 
inductors  located  between  the  elements  are  designed  to  saturate  at  the  conclusion  of  the 
charging  pulse  thereby  connecting  all  the  elements  of  the  line  together  and  allowing  the 
RF  wave  to  propagate  down  the  line  to  the  output.  This  frozen  wave  generator  design  is  a 
bit  more  complicated  than  the  version  described  below  but  does  allow  for  cen te”  charging 
and  eliminates  transformer  coupling  problems. 
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The  second  concept  illustrated  in  Fig.  5.15  uses  saturable  transformers  between  ele¬ 
ments.  The  transformers  are  used  to  create  the  wave  on  the  line  while  in  the  unsaturated 
state.  At  the  conclusion  of  the  charging  pulse,  the  core  material  making  up  the  trans¬ 
formers  is  designed  to  saturate.  Once  this  saturation  occurs,  the  line  is  decoupled  from 
the  nonlinear  magnetic  charging  system  and  the  saturated  transformers  now  act  els  the 
inductors  of  a  lumped  element  transmission  line. 
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